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1. Introduction 

Parton shower Monte Carlo event generators, such as Herwig [1] and Pythia [2], have 
proven to be enormously useful since the development of the main ideas in the 1980s 
[3, 4, 5]. These computer programs perform calculations of cross sections according to an 
approximation to the standard model or some of its possible extensions. Because of the 
great success of these programs, it is worthwhile to investigate possible improvements. In 
this paper, we propose a theoretical structure for event generators that generalizes the 
structure of current programs and allows the elimination of certain approximations used 
currently. 

Parton showers are mostly reflections of QCD interactions. In order to present a 
reasonably complete discussion of the QCD issues in a parton shower while keeping the 
length of this paper within reasonable bounds, we limit the presentation to QCD and omit 
any discussion of how electroweak and beyond-the-standard-model interactions are to be 
added to the QCD interactions to make a useful event generator. 

What is a parton shower Monte Carlo event generator? Let us consider hadron-hadron 
collisions, which is the case relevant for the Tevatron and the Large Hadron Collider. An 
experiment will produce a large number of events /, where one can characterize an event as 
a list of the momenta and flavors of the final state particles produced. The experiment can 
measure a cross section (t[F] corresponding to an observable"^ that assigns to each event / 
a number F{f). The relation of the cross section and the function F is 



where C is the integrated luminosity for an experimental run and the fn are the observed 
events. For example, the cross section to produce a Higgs boson and two jets having certain 

^In order to be subject to reliable calculation in QCD perturbation theory, the function F should have the 
property known as infrared safety. However, a parton shower event generator is also useful for observables 
that are not infrared safe. 
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characteristics is specified by setting F(f) = 1 if / contains a Higgs boson and two jets 
having these characteristics and F{f) = otherwise.^ 

A parton shower Monte Carlo event generator calculates this cross section by producing 
a large number N of simulated events fn, each with an accompanying weight Wn- The 
calculated cross section is then 



Most typically, the weights are all equal, so that l/w^ is the simulated luminosity per point 
C/N. Our definition of the category of parton shower Monte Carlo event generator includes 
the possibility that the weights are complex numbers produced for each event. It is always 
possible to throw away the imaginary parts of the Wn since we know in advance that the 
imaginary part of the sum in Eq. (1.2) vanishes, so having complex weights is equivalent 
to having real weights that can be positive or negative. This situation occurs in typical 
event generators [6, 7] that are based on next-to-leading order perturbation theory.'^ 

In a typical parton shower event generator, the physics is modeled as a process in 
classical statistical mechanics. Some number of partons are produced in a hard interac- 
tion. Then each parton has a chance to split into two partons, with the probability to 
split determined from an approximation to the theory. Parton splitting continues in this 
probabilistic style until a complete parton shower has developed. 

The parton splitting probability is biggest when the two daughter partons are almost 
massless with nearly coUinear momenta or when one of their momenta is soft (near p = 0), 
or both. There is a simple underlying approximation used: the amplitude for producing 
m + 1 partons when two of the momenta pi and pj are nearly collinear or one is soft factors 
into a splitting function times the matrix element for producing m partons. 

The underlying approximation is the factorization of amplitudes in the soft or collinear 
limits. However, further approximations are usually added: 

1. The interference between a diagram in which a soft gluon is emitted from one hard 
parton and a diagram in which the same soft gluon is emitted from another hard 
parton is treated in an approximate way, with the "angular ordering" approximation. 

2. Color is treated in an approximate way, valid when — > where A'c = 3 is the 
number of colors. 

3. Parton spin is treated in an approximate way. According to the full quantum ampli- 
tudes, when a parton splits, the angular distribution of the daughter partons depends 
on the mother parton spin and even on the interference between different mother- 
parton spin states. This dependence is typically ignored. 

^The case in which F{f) takes values or 1 is the most common, but other possibilities are allowed. 
For instance, the energy-energy correlation function in electron-positron annihilation is of the more general 
variety. 

^The recent paper [8] provides an exception to this rule. 





n=l 



- 3 - 



With the use of these further approximations, one can get to a formahsm in which the 
shower develops according to classical statistical mechanics with a certain evolution oper- 
ator. 

Our purpose in this paper is to investigate whether one can have a formulation of 
parton showers based on the factorization of amplitudes in the soft or collinear limits in 
which one does not make the additional approximations enumerated above. For this, one 
would have to use quantum statistical mechanics instead of classical statistical mechanics. 

It might seem that doing the problem in quantum mechanics is hopelessly complicated. 
However, within the soft /collinear factorization approximation, the problem is fairly simple 
because it is almost classical. In fact, if partons did not have color or spin, the problem 
would be classical (as we discuss in Sec. 2). Thus what we need is a fully quantum treatment 
of color and spin. We arrange for this by making use of the quantum density operator in 
color (g) spin space. 

In the subsequent sections, we define evolution equations for the quantum density ma- 
trix within the soft /collinear factorization approximation. The matrix evolves in "shower 
time" from harder splittings to softer splittings. The iterative solution of these equations 
gives cr[F] in the form of a sum of integrals. To give some idea of the structure, we omit 
any mention of hadronization and write the result in a notation that is quite abbreviated 
compared to the notation in the body of the paper. 




There is, first of all, an integration (including sums, for discrete variables) over momenta, 
flavors, spins, and colors for initial partons that emerge from the hard matrix element 
and its complex conjugate. Here we call all of these variables collectively Pq, the initial 
partonic variables. There is a function /o that depends on Pq and represents the hard 
matrix element at the start of the shower times its complex conjugate. Then there is a 
sum over how many splittings, N, there are.^ Next there is an integration over splitting 
variables Cj for the jth splitting. The splitting variables include the label telling which 
parton split and momentum variables, for which a dimensionless virtuality y, a momentum 
fraction z, and an azimuthal angle (j) might be used. There are also discrete flavor, color, 
and spin variables. At each splitting, there is a set of starting partonic variables, Pj-i 
and a set of new partonic variables Pj that are determined by Pj-i and the splitting 
parameters Cj- For each splitting, there is a function fj that depends on Pj-i and (j. We 
have integrations over the splitting parameters for splittings 1 through A^. At the end, 
there is the measurement function F that depends on the partonic variables P/v reached 
after all of the splittings. 

The structure of this representation is similar to that in conventional parton showers, 
with the functions fj made from splitting functions and Sudakov exponentials that express 

*We have formally iterated the evolution equation an infinite number of times, allowing any number of 
splittings. However, we imagine that there is a cutoff on splitting hardness, so that very large values of A'^ 
are seldom encountered. Some of our splittings are 1 — > 1 self interactions rather than 1^2 splittings. 
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the probability for not splitting. There are, however, some important structural differences 
that result from not making the approximations 1, 2, and 3 above. Chief among them is 
the use of the spin and color variables. 

What we develop in this paper is an evolution equation that results in a representation 
of (t[F] as integrals of known functions. Of course, one will want to turn the integrals into 
numbers. How to do that is a question of numerical integration that we leave for future 
work. However, it may be useful to sketch how a numerical evaluation might work. 

To evaluate (t[F] numerically, one has to construct the functions involved as described 
in the body of this paper and then perform all of the integrations (and sums). In particular, 
Monte Carlo integration can be used for many of the integrations. In a numerical method 
that is very, very simple, one can choose random points Pq first, according to a density 
po- Then one would choose the first splitting variables, Ci, according to a density pi 
that is determined by the variables Pq. This would determine new partonic variables Pi. 
Continuing, we choose C,j according to a density pj that is determined by variables Pj-i 
and we use Qj together with Pj-i to determine Pj. At the end, we use P/v as input to the 
measurement function F. This process constitutes a Markov chain that produces "events" 
with a final state P/v- The probability density for getting a final state P/v with a given 
shower history is /?o times the product of the pj for j > 1. We multiply F{Pn) by a weight 
equal to /0//O0 times a product of the fj/pj- 

It would be a design goal to choose the pj to be roughly proportional to the absolute 
values of fj. This kind of importance sampling would produce weights that do not vary 
over a wide range. 

We may note that in a conventional parton shower, the fj are everywhere positive and, 
for J > 1, integrate to 1. Thus one can choose the po to be proportional to /o and pj = fj 
for J > 1. Then the weight function is a constant. In our case, the factors in fj are not 
everywhere positive, so we expect to need weights, which could have either sign. 

What we have described above would generate a very conceptually simple numerical 
solution to the evolution equation. We expect that one could do much better, particularly 
by performing the spin sums not by numerical Monte Carlo summation but by exact 
summation. For this, one could adapt the method proposed by Collins [9] and elaborated 
by Knowles [10] and by Richardson [11]. 

We leave issues of the numerical evaluation of the integrals, beyond this simple discus- 
sion, for future work. 

Some features of the formalism presented here can best be understood by asking what 
would happen if we kept only the leading 1 /N'^ — > limit and averaged over spins every- 
where, thus making approximations 2 and 3 above. We would then have a shower based on 
gluon emission from color dipoles. With such a picture, the imposition of a cut to enforce 
angular ordering (approximation 1) is not needed: interference between gluon emissions 
from both halves of a color dipole is already included. The resulting evolution equation 
would then be similar to what is implemented as A;_L-showers in Pythia, as described in 
Ref. [12]. One could also modify the formalism presented here change from k"^ as the evolu- 
tion variable to the version of k\ used in Pythia. Then there would be two main features 
that differed between the present formalism and Pythia. One difference is in the choice 
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of splitting functions. The Pythia choice is the Altarehi-Parisi sphtting functions, defined 
with a certain definition of the momentum fraction z. Our sphtting functions are made 
from the Feynman diagrams for one off-sheh parton producing two on-shell partons with 
only a minimal manipulation to separate this part of the diagram from the hard scattering 
to which it attaches. The definitions match in limit of collinear splittings, but differ away 
from this limit. The other difference is in the momentum mapping that connects the mo- 
mentum space for two initial state partons and m + 1 final state partons to that with m 
final state partons. We have investigated [13] the possibility of using the Catani-Seymour 
[14] mapping and splitting functions, which are commonly used for next-to-leading order 
calculations. However, we have here adopted a mapping that avoids the use of designated 
"spectator" partons that share some of their momenta.^ Pythia uses the Catani-Seymour 
momentum mapping for final state splittings and something more complicated for initial 
state splittings. We note here the recent paper [16], which explores other possibilities for 
both momentum mapping and splitting functions. We also note that S. Schumann and 
F. Krauss and, separately, M. Dinsdale, M. Temick, and S. Weinzierl have very recently 
implemented the Catani-Seymour dipole subtraction functions and momentum mappings 
as the basis for a parton shower. The first results look promising [17]. 

One of our goals has been to have a formulation that can coexist easily with matching 
the probabilities generated by showering to known exact tree level matrix elements, as in 
Refs. [18] and with using next-to-leading order hard matrix elements as in Refs. [6, 7, 19]. 
However, we leave for future work the analysis of how one can match the showers to the 
exact tree level matrix elements or to next-to-leading order hard matrix elements. 

One may wonder whether removing the approximations 1, 2, and 3 listed above is 
numerically important. We do not have a definitive answer. What we would like to do is 
to set up a formalism that does not make approximations beyond the basic soft /collinear 
factorization approximation, then (in future work) implement this formalism as a working 
algorithm and computer code. One could then make the further approximations separately 
or all together and see what difference they make. 

We close this introduction with some comments on whether a parton shower Monte 
Carlo event generator ought to allow weights for generated events, and in particular nega- 
tive weights. 

We first note that in a real experiment the relation between the measurement function 
F and the measured cross section is a little more complicated than we indicated in Eq. (1.1). 
Instead, we have 



where £ is the integrated luminosity for an experimental run in which N total events n 
are collected and a.„ is the acceptance for the event resulting from the way the detector 
is triggered. For example, if a fraction 10^'^ of a certain kind of event is recorded, then 

^One motivation for avoiding a special role for designated "spectator" partons is that if one wants to go 
to a next-to-leading order splitting kernel or to subtractions for a NNLO perturbative calculation, problems 
can arise from a third parton becoming collinear with the designated spectator parton [15]. 
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for those events is 10 ^. Thus the weight factors Wn/N in Eq. (1.2) are analogous to 
l/{anC) in the analysis of real data. 

We can also examine the effect of weights on the statistical error in Eq. (1.2). The 
expected error £ is given by 



Here (wF) indicates an average. If it is very expensive to use a large N, for instance 
because calculating F{fn) requires a full detector simulation, then one would like to make 
S"^ as small as possible for a fixed A^. That suggests not that Wn should be constant, but 
that Wn F{fn) should be approximately constant for the observables F of most interest.^ 
It is never possible to make Wn F{fn)/ {wF) = 1 for all events, but one does not want to 
have lots of events for which this ratio is much smaller than 1 nor any events for which the 
ratio is much bigger than 1. Having events for which Wn F(fn)/ (wF) ~ — 1 is not a good 
thing, but, since | — 1| is not much larger than 1, it is not really damaging from the point 
of view of avoiding large statistical errors. 

Evidently having weights that are real numbers of either sign, or complex numbers, 
does not make it impossible to apply Eq. (1.2). It does, however, make the analysis a little 
more complicated. However, we believe that that the added complication does not present 
a serious problem. 

Preview. Since this is a rather lengthy paper, some preview of what is in it may be 
helpful. Section 2 contains an introduction to the notation we use. This notation is, 
we think, useful for thinking about a variety of formulations of the parton shower idea. 
We present it in the context of a simple scalar field theory that is free from a lot of 
the complications of quantum chromo dynamics (QCD). We then turn to QCD, with its 
complications. We present in Sec. 3 the structure that we propose for a parton shower that 
contains quantum interference. In order to present this structure in just a few pages, we 
leave for later sections most of the detailed definitions. The first of these, the momentum 
and flavor mapping, is covered in Sec. 4. Then Sec. 5 covers spin. This provides enough 
background to present the splitting functions for the quantum amplitudes in Sec. 6. The 
description of color, which is rather more complicated than that of spin, is presented in 
Sec. 7. We are then able to specify the shower evolution operator in Sees. 8 and 9. We 
follow this with discussions of two interesting issues, the evolution of color in Sec. 10 and 
soft gluon coherence in Sec. 11. We analyze the structure of the functions that appear in 
the Sudakov exponent in Sec. 12. Eventually shower evolution stops and a hadronization 
model is inserted. We discuss this in Sec. 13. We present some concluding remarks in 
Sec. 14. There are two appendices that deal with certain technical issues. 

®Thus, if our primary interest were in the high Pt tail of a jet Pt distribution, we would not want to 
use most of the available computer time to generate low Pt events. Rather, we would want to generate few 
low Pt events, giving each of them a high weight to compensate. 
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2. A notation for pctrton showers 

Starting in the next section, we present a formulation for parton showers in QCD hard 
scattering events in hadron-hadron collisions, taking into account the complexities intro- 
duced by spin and color correlations and by soft, wide angle gluon emissions in addition 
to collinear spUttings. In order to do tliis, we use a mathematical language that helps 
to organize the algorithm. Alas, the complexities of the real physical situation make the 
needed construction a bit subtle. Therefore, in this section we first introduce some of the 
needed language in a simpler situation. The notation introduced here will be used again 
for QCD in the subsequent sections. 

Consider the process e~^e~ hadrons in a world in which hadrons consist of just one 
kind of massless scalar particle, which has no color. The e^e^ annihilation produces a pair 
of virtual scalar particles through an interaction that we need not specify. In the evolution 
of the hadronic state, we can still have collinear singularities similar to those found in 
QCD if the theory consists of (f)^ theory in six dimensional space-time. The cross section 
to measure an observable F can be written as 

^[^1 = Y.hl K^'}-] \M{{p}„^)\'F{{pU) . (2.1) 

TTl. J 

m 

Here there is a sum over the number m of produced particles, {p\m = {pi, ■ • • ,Pm}, and 
j [d{pU] ^ n {/ ^ 27r<5+(pf )| {2t,)H(p, -Y^p^ , (2.2) 

with Po = iVs, 0). The function M{{p}m) gives the matrix element to produce m particles 
with momenta {p}rn, while F({p}rn) describes the measurement to be done. 

We want to describe this using an algorithm that approximates |M({|)}„j)p based on 
|M({p}^)P for m = 2 and the subsequent generation of the rest of the particles based on 
a narrow angle approximation for one particle to split into two. This is to be done using 
a Monte Carlo simulation in which the system evolves from 2 particles to many particles 
as a simulation time t progresses from t = to a large value, at which the simulation is 
terminated. There are various possibilities for the physical meaning of the time t. We will 
take it that a splitting I ^ i + j occurs at time t = \og[Q'^/ {2pi ■ Pj)), where Qq is the 
hardness scale of the hard interaction with which we start. 

At each stage of this simulation, let the cross section to have m particles with momenta 
{p}m be p{{p}m,t). Summing over the number of particles and integrating over momenta 
gives the total cross section, 

^^ = Y.—\ I [^Mm] P{{p}m,t) . (2.3) 
TTl. J 

m 

At the final time, , the value of the measurement function is 

'^[^l = J2^J id{p}m] Pi{p}m, if) FHpU) . (2.4) 
m ' 
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The possible functions p (at a given time t) form a vector space, so that p at time 
t can be considered to be a vector \p{t)). We use rounded brackets here. The notation 
\'4)) is reserved for a quantum state, while |p) denotes a state in the sense of statistical 
mechanics. We therefore call it a statistical state. The inner product is^ 

(^1^) =T.hl ^(W-) ^(W-) • (2-5) 

m 

We can define basis vectors | {p}m) in this space so that 

P{{p}m,t) = {{p}m\p{t)) . (2.6) 

With these definitions, there is a completeness relation 

^ = J2^ [ \{p}m) {{p}m\ • (2.7) 

III. J 

m 

The measurement function F can also be considered to be a vector, Thus 

a[F] = {F\p{ti)) . (2.8) 

There is a special vector (l| with 

(1|M„) = 1 . (2.9) 

This vector represents the totally inclusive measurement function corresponding to the 
total cross section, 

(7T=(l|p(i)) • (2.10) 

Now we are ready to discuss the evolution of the statistical state. We take the evolution 
to be given by a linear operator U{t,t'), with 

\p{t))=U{t,t')\p{t')) . (2.11) 

Here U{t,t) = 1. These operators have the group composition property 

U{h,t2)U{t2M)=U{hM) ■ (2.12) 

The class of evolution operators that we will use is defined by two operators. The first 
is an infinitesimal generator of evolution or hamiltonian, Hi{t). We can specify 'Hi{t) by 
giving its action on an arbitrary state \p), 

niit)\p) . (2.13) 

In a lowest order shower, the operator 7ii{t) describes parton sphtting, changing a state 
with m particles to one with m+ 1 particles. One of the particles in {p}m, say particle I, is 



Note that there is no * here. 
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removed and replaced by two, with momenta pi and Pm+i- The quantum amplitude after 
the splitting is approximately 

M({p}„+i) « M{{p}m) X —4 . (2.14) 

^Pl -Pm+l 

Here, to precisely define the right-hand side, one needs to redefine the momenta so that the 
mother parton has a momentum pi ~ pi + p^+i that is nevertheless on-shell, p^ = 0. Thus 
the momenta {p}m are functions of the momenta {p}rn+i- We omit a discussion here of the 
various ways to define this momentum mapping. More important for now is the idea that 
the quantum matrix element factorizes in the form (2.14) when 2pi ■ Pm+i is much smaller 
than any of the dot products among the momenta in {p}m- This factorization is at the 
heart of the reason why parton shower Monte Carlo programs give useful approximations. 
For the statistical splitting function in Eq. (2.13), we need the square of the quantum 
amplitude. Thus we want 

2 

{{p}m\p) ■ (2.15) 

We have inserted the definition of the Monte Carlo time, t, that we here imagine using. For 
our present pedagogical purposes, the details of the definition of 7Yi(i) are not so important. 
What is important is that it reflects the factorization (2.14). 

The second operator used in the construction oiU{t' ,t) is a no-change operator M{t, t') 
with Af{t, t) = 1 and 

AA(t3,i2)AA(i2,ti) =AA(t3,ti) . (2.16) 

The no-change operator leaves the basis states unchanged except for multiplying each of 
them by an eigenvalue A: 

M{t',t)\{pU) = A{t',t;{pU)\{pU) . (2.17) 




Figure 1: Illustration of Eq. (2.18). The wide ovals represent the evolution operator U while the 
narrow ovals represent the no change operator A/", which provides the Sudakov exponentials, and 
the small circle is a parton splitting. Shower time t runs from left to right. 



The evolution operator lA is expressed in terms of the hamiltonian and the no-change 
operators by (see Fig. 1) 

Uit3,ti) = M{t3,h) + dt2U{h,t2)'Hl{t2)M{t2M) ■ (2.18) 

This equation is interpreted as saying that either the system evolves without splitting from 
tl to ts, or else it evolves without splitting until an intermediate time t2, splits at t2, and 
then evolves (possibly with further splittings) from t2 to t^. 

We need one more ingredient. We wish to construct the evolution so that it leaves the 
total cross section invariant, {l\U{t' ,t)\p) = (l|/o). This should hold for every state j/o), so 

{l\U{t' ,t) = {l\ . (2.19) 

This assumption defines what AA(t3,ti) has to be. Multiply Eq. (2.18) on the left by (l| 
and on the right by |{p}m)- Then, using Eqs. (2.19) and (2.17) we have 

l = A(t3,ti;Mm)+ / dt2 (l|Hl(t2)|Mm) A(t2,ti;Mm) . (2.20) 

If we now differentiate with respect to ^3, we get 

-^A(t3,ti;Mm) = -(l|Wl(t3)|Mm) A(t3,ti;Mm) ■ (2.21) 

The solution of this with the initial condition A(ti,ti; {p}m) = 1 is 

A(t3,ti;Mm) =exp(^-^'dr (l|Hi(r)|{p}™)^ . (2.22) 

This result has a simple interpretation. The quantity (l|'Hi(r)|{p}m) is the total proba- 
bility for one of the partons in the state {p}m to split at time r. The exponential, known 
as the Sudakov factor, is the probability that none of these partons splits between ti and 
ts- 

This operator notation provides a convenient way to express the essence of standard 
shower Monte Carlo algorithms. Different algorithms differ in how the evolution variable t 
is defined and in exactly what the splitting operator 7ii{t) is. 

At its heart, the shower Monte Carlo idea is that soft interactions factor from hard 
interactions in QCD. After being produced in a hard interaction, partons in QCD travel 
a long way before undergoing much softer interactions such as splitting. What subsequent 
splitting does occur does not much change the basic hard matrix element for producing 
(almost) on-shell partons. There are some complications, however. On shell partons carry 
both spin and color. Furthermore, soft gluons can transmit color changes over long dis- 
tances. In the subsequent sections, we will extend the meaning of the symbols used here 
so as to accommodate spin and color. 

3. Structure of the calculation 

In this section we introduce notational conventions and a general structure for the calcu- 
lation that we will use later in the paper. 
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3.1 The space of quantum parton states 

In order to describe showers, we need a notation for the description of quantum states 
consisting of two initial state partons and m final state partons. The partons are labeled by 
an index such as i that takes values "a" or "b" for the initial state partons and 1,2, ... ,m 
for the final state partons. Each parton is described by a momentum p, a flavor / G 
{g,u,u,d, d, ...}, a spin index s and a color index c. We denote the quantum numbers of 
such a state by 

{P,f,S,c}m = {??a,a,Sa,Ca;r?b,ft, Sb,Cb;Pl,/l,Sl,Ci; ...■,Pm,fm,Sm,Cm} ■ (3.1) 

Here there is a special notation with respect to the incoming partons. The momentum 
fractions of the incoming partons are denoted by r/a and rj]^, defined below. We call the 
flavor of parton "a" simply a and we call the flavor of parton "b" simply b. This notation 
is useful for designating the parton distribution functions, fa/AiVa.^ I^f) ^^"^ fb/BiVh, fJ'p)- 
For the purpose of describing backward evolution of the initial state partons, we will often 
need the antiflavors of the incoming partons. We use /a and /b for these, 

= - " ' (3.2) 
h = - b . 

Here our notation is — u = u, — u = u, — g = g, etc. 

The final state partons are always on-shell. Our kinematics allows parton masses, with 

p]=m\fj) . (3.3) 

To describe the momenta of the initial state partons, we start by defining pA and pB to be 
massless approximations to the momenta of the two incoming hadrons 

pa = o , 

pI = , (3.4) 
2pA-PB = S . 

The initial state partons are on-shell. In general, they can have masses but, with a small 
modification of the notation, their masses could be set to zero.^ In any case, we take the 
initial state partons to have zero transverse momentum. Thus 

Pa. = VaPA H PB , 

VaS 

, m^(/b) 
Ph = VhPB H PA ■ 



(3.5) 



This defines the momentum fractions r^a and Tjy,. 

In the event that we include parton masses, the hardness scale Qq at which the parton 
shower is initiated should be much bigger than the mass of any parton that is included 



*There are several possibilities for the treatment of masses of initial state partons and there are some 
subtle issues associated with the choice. We mention some of these issues in Sec. 14. 
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as a possible constituent of the incoming hadrons. For example, if the hard process were 
t + t production near threshold, then top quarks should not be used as possible initial state 
partons. Thus we demand that 

Ql > ^rnl , (3.6) 

where mn is the mass of the heaviest quark that is included as an initial state parton, 
typically the b quark. In any reasonable application of the formalism of this paper, the 
">" here will be "S>." Equation (3.6) suffices to make certain kinematic formulas in the 
paper work. At any stage in the shower, we define Pa + Ph = Q- At the first step, the 
starting hard scattering, with a sensible definition of the starting hardness scale we must 
have > Qq. The values of increase as the shower develops. Thus at any stage we 
will have 

(Pa + Pb)^ >4m| , (3.7) 
We will impose a kinematic restriction on the momentum fractions. 

These limits require that the momentum in the system in the direction of pA comes mainly 
from parton "a," and the same for pA ^ PB, a ^ b. We can be sure that both conditions 
hold by requiring 

??a??bS > m\ . (3.9) 

The reason for imposing this condition is as follows. Given Eq. (3.5), there are two choices 
for r]ar]hS that yield the same value of (pa + Ph)'^- As long as condition (3.7) holds, the 
larger of the two choices for riarj\^s satisfies Eq. (3.9). Imposing Eq. (3.9) eliminates the 
other solution, in which parton "a" moves in approximately the p-Q direction and parton 
"b" moves in approximately the p\ direction, creating a large value for {p^ + Phf'- 

The upper limit on "q^ and % is 1. We note here that this is an approximation. To 
discuss this, let nA and n-Q be dimensionless lightlike vectors in the directions of pA and 
Pb respectively, normalized to ua ■ n-Q = 1. Then consider, for example, the limit on r/a. 
The total momentum in the direction of riA of the final state particles is 



(Pa +Pb) • "-B 



PA-riB ■ (3.10) 



This can be bigger than the available momentum pA ■ ^^B if Va. is very close to 1. A remedy 
for this would be to redefine pA and pB in our formulas. Suppose that the exact hadron 
momenta are Pa and Pb- Then if we put 

PA = A Pa • "B ^A , 
PB = X Pb ■ nA nB , 

then there is a value of A that makes (pA + Pb)^ = (Pa + Pb)^- By taking a value of A 
that is a little smaller than this, one can ensure that the momentum in the final state in 
the directions of ua and hb is not more than was present in the initial state. It is this 
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momentum that is available for the "underlying event." With this adjustment, the value 
of s in our formulas is a little less than the true cm. squared energy, (Pa + Pb)^- 

Our notation with respect to spin and color is meant to be flexible. A standard helicity 
basis will work for spin. For color, we begin with a straightforward basis in which each 
parton i has a color index Cj that can take values 1,2,3 for quarks and antiquarks and 
1,. . . ,8 for gluons. Later, we will want to consider the subspace of the whole color space in 
which the parton state is a singlet under the SU (3) color group. Wc will choose a basis for 
this subspace. Using this basis, we will still have labels that we can call {c}m, but the new 
labels will describe the color links among the partons rather than individual color indices 
Cj for the individual partons. The notation {p, f, c, s}m is thus supposed to include the 
possibility of any representation of the colors of the m + 2 partons. 

3.2 The density matrix 

A matrix element used in the computation of a cross section can be thought of as having 
the form of a function of the momenta and flavors that carries indices for spin and color, 

M({p,/},„)^:;^-,?i::.fc • (3-12) 

Here we denote the functions M for different numbers m of final state partons by the same 
name, simply M({p, /}m) rather than Mjn{{p, f}m)- The array M can be thought of as a 
vector in spin and color space, 

\M{{pJU)) . (3.13) 

The inner product (M'|M) denotes multiplying M'* by M and summing over the spins 
and colors. 

An observable F can be specified by giving a set of functions F{{p, f}m) that are linear 
operators on the color-spin space. (In many important cases, F{{p, f}m) is a function times 
the unit operator on color-spin space. If F is be simply made from theta functions defining 

final state cuts, then (t[F] is the cross section to find the final state partons within the 
cuts.) With this notation the cross section for an observable F takes the form^ 



f a/ A{r)a., IJ-l) fb/ siVh, 14) 

4nc(a)nc(6)2?7ar/bPA-PB (3.14) 



m 

X (M({p, f}m)\F{{p, f}m)\M{{p, f}m)) 



Here the functions / are parton distribution functions while nc(a) is the number of colors 
that a parton of fiavor a can have, iVc = 3 for a quark or antiquark, A^^ — 1 = 8 for a gluon. 
The factor Anc{a)nc{b) turns the sum over spins and colors for the initial state partons 
into an average over spins and colors. We have indicated the appropriate integrations over 



^This formula contains a parton flux factor 2ris,rihPA ■ Pb that corresponds to massless partons. The flux 
factor for scattering of free massive particles is more complicated. However, (pa+Pb)^ is always bigger than 

the initial hard scale Qq and the formalism of this paper is valid only when is much larger than the 
masses of any initial state partons. For this reason, we use the flux factor for massless parton scattering. 
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momenta by 



X (27r)^5|^Pa +Pb - E^O ^("^H < ^a%s) 
X5({P,/}m) • 



(3.15) 



Here g{{p, f}m) is an arbitrary function. 

The final state particles carry labels i £ {1, • • . Then particle i has momentum, 

flavor, spin, and color given by {pi, fi,Si,Ci}. One can arrange the definitions such that the 
amplitude M is symmetric under interchange of the labels. However, we do not necessarily 
do so. Instead, the notation allows for a general labeling scheme. A measurement function 
F must be symmetric under interchange of labels, since the labels are not physical. Two 
amplitudes M that become the same if the labels are symmetrized are equivalent. 

We will find it useful to rewrite cr[F] in the form of a trace over the spin and color 
space, 

^[^] = E^ / [d{pJU]TT{pi{pJ}m)F{{pJU)} , (3.16) 



where 



(3.17) 



Thus p is the density operator in color spin space. It is illustrated in Fig. 2. The density 
operator, for momentum as well as spin, is widely used as the basis of quantum statistical 
mechanics. It was introduced for the spin space in parton showers by Collins [9] and is 
used in Herwig for the heavy partner particles in supersymmetry [11]. 

We can expand p{{p, f}m) in basis states |{s,c}„i) for the color (g) spin space, 

P({P, /l^n) = EE K^''^^™)^('LPi/>5'>c',S,c}m,) ({s',c'}m,| ■ (3.18) 
s,c s',c' 

Here p{{p, /, s' , c' , s, c}m) is a function that depends on the momenta and flavors {p, f}m, 
the labels {s,c}m for the quantum "ket" state and the labels {s',c'}m for the conjugate 
quantum "bra" state. We denote the state labels collectively by {p, f, s' ,c' , s,c}m- Thus 
p{{p, f, s' , c' , s, c}m) a function giving the matrix elements of the density matrix. We find 
it convenient to base our treatment on this function. 



^''Just to take a trivial example, in a utig state, the label 1 might be assigned to the up quark, 2 to the 
anti-up quark and 3 to the gluon. Of course, it will not work in general to use the flavors as labels because 
one can have two final state partons with the same label. 
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Figure 2: Illustration of a contribution to p, Eq. (3.17). The Feynman graph on the left is a 
contribution to \ AI {{p, ,f}m)) and the Feynman graph on the right is a contribution to (^l\I{{p, f}m) \ ■ 
The factor with parton distributions is not illustrated. The momenta and flavors of the labeled 
partons match between |M({p, /},«)) (A/({p, /},„) | , but the partons on the left have colors and 
spins {s, c}m while those on the right have possibly different colors and spins {s', c'}m- Additionally, 
p can contain quantum interference between different graphs, as illustrated here. We imagine that 
eventually the partons will evolve to form a final state in the middle and an initial state on the left 
and right. 

Later, we will choose an orthonormal basis for the spin space, so that 

{{s'}m\{sU) = 6l'^l- . (3.19) 

For the color space, we will find it convenient to use a basis in which ({c}m|{c}m) is 
approximately but not exactly 1 and where ({c'}m|{c}m.) is not generally zero for {c}m. / 
{c'}m- With a non-orthogonal basis, we need to be a little careful about the notation. 
Suppose that we expand a vector in such a basis, 

= E \Mm)a{{c}m) . (3.20) 

{c},„ 

A convenient representation for the expansion coefficients a({c}m) is obtained by taking 
matrix elements with elements of the dual basis |{c}.m)^ defined^^ by 

^({c'}^|{c}„) = 4Jj . (3.21) 

Using the dual basis we can write 

^({c}m|^) = a({c}„^) . (3.22) 

The expansion of any vector in the basis |{c}m) can conveniently be obtained directly by 
using the completeness relation 

1 = E K^}-) d({'^}™I • (3.23) 

{c}m 

^^Readers familiar with general relativity can think of ({c}mjci) as the covariant components flc of a 
and ^(^{c}m\ci) as the contravariant components a'^. One can lower indices using the metric tensor g^c' , 

analogous to {{c}m\{c'}m)- The inverse matrix {{c}m\{c'}m.) is analogous to g'^'^ . 
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If we want to expand a bra vector in the basis ({c}m| we can use the completeness relation 

1 = E K^}->D (i'^}™! • (3.24) 

This is particularly useful with respect to operators. Let O be a linear operator on the color 
space. If 1^) has expansion coefficients a{{c}m) defined by Eq. (3.20) and = 0|^/^) has 
expansion coefficients a'{{c'}m) then O is conventionally described by the matrix defined 

by 

a'{{c'}m) = o{{c'}m,{c}m)a{{c}m) ■ (3.25) 

{c}m 

We can write this as 

0{{cU,{cU) = J,{c}m\0\{c}m) ■ (3.26) 

3.3 Statistical states 

The set of functions p{{p, f, s', c', s, c}m) forms a vector space, which we can call the statis- 
tical state space (as distinct from the quantum state space). We can call the vector corre- 
sponding to this function simply | /?) . Note the rounded brackets instead of angle brackets 
that we use for quantum states, as in |^). We can define basis vectors ({p, /, s', c', s, c}m\ 
for the statistical state space-^^ so that 

P{{P, f, s', c', S, C}m) = {{P, f, s', C, S, C}m\p) ■ (3.27) 

There are also ket basis vectors such that the completeness relation for the basis states is 

1 = E;^/ [d{P^f^^'^(^'^^^^}m]\{P,f,s',C,S,c}m){{P,f,s',C,S,c}m\ , (3.28) 

where [d{p, f, s, s', c, c'}m\ is an extension of the integration measure Eq. (3.15) 

[d{pJ,s',c',s,cU]^ J [d{p,f}m] E E n| E f • (3-29) 
The corresponding inner product of basis states is 

{{pJ,s',C,S,c}rn\{p,f,S,C,S,c}rh) = Sm,rh 6{{p, f , s' , C , S , ajm] {p, f , S , C , S , c}m) , 

(3.30) 

where the function 5 is a generalization of the ordinary delta-function that is defined by 

I [d{p, f,S ,C ,S,c}m] 5[{pJ,s\c ,S,c}m\{pJ,S ,C ,S,c}m) h{{p, f, s' , C , S, c}m) 
= H{P,f,S,C,S,c}m) ■ 

(3.31) 



^^More precisely, the bra vectors (Fj are vectors in the dual space to the ket vectors |p), that is the space 
of linear functions on the vectors jp). 
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Here h is any well behaved function of the variables indicated, defined on the integration 
surface. 

Let us define a vector corresponding to a measurement function F using^'^ 

{F\{p,f,S,C,S,c}m) = {{s',c}rn\F{{p,f}m)\{s,c}m) ■ (3.32) 

Then, using the completeness relation (3.28), the cross section (3.16) corresponding to a 
measurement function F can be expressed using Eqs. (3.18) and (3.27) as 

a[F] = {F\p) . (3.33) 

In the case that F consists of a function F{{p, f}m) times a unit operator in spin-color 
space, the inner product {F\p^ is 

i^\p) = y]~7 / [d{P,f,s',c',S,c}m]F{{pJ}„^) {{s'}„i\{s}m){{c}m\{c},n) 

' ^ ' ^mlj ' (3.34) 

X P{{pJ,s',c',S,c}m) . 

3.4 The resolution scale 

We now need to introduce a resolution scale into our equations. We first discuss the reso- 
lution scale of the observable. In Eq. (3.33), let us suppose that the observable represented 
by the functions F{{p, f}m) is infrared safe. To be precise about what this means, we first 
demand that the functions F{{p, f}m) be smooth functions of the momentum variables and 
that they be invariant under label interchanges. Then we consider a list of parton variables 
{p, f}m+i for m+1 partons and suppose that Pm+i becomes collinear with the momentum 
Pi of parton I. Then we can consider the list of parton momenta {p, f}m where pj = pj 
and fj = fj for j ^ I while pi = pm + Pi and fi = fm + fi (with the obvious definition 
of adding flavors). That is, the partons with variables {p, f}m+i could have arisen from 
partons with variables {p, f}m by the collinear splitting of parton / into new partons with 
labels / and m + 1. Then infrared safety requires that 

Hip, /}m+i) ^ F{{p, fU) (3.35) 

in the limit. This should also hold if parton m + 1 becomes soft, Pm+i 0. It should also 
hold with a suitable adjustment of the notation for a splitting of one of the initial state 
partons. This is, so far, just the standard definition of infrared safety. It allows us to have 
perturbatively calculable cross sections. Now let us extend the definition to include a scale. 
We can say that the observable is infrared safe at scale fi^ if 

F{{pJU+i)^F{{pJU) (3.36) 

when \2pm+i -pil < /U^. To be really precise, we should specify how good this approximation 
has to be, but this will not matter for our purposes. Strictly speaking an "infrared safe 

^''Note that this equation for the measurement function has a different structure from the equation used 
to define the statistical state vector |p), p({p,/}m) = \{s, c}m){{p, f, s' , c' , s, c}m\p) {{s' ,c'}m\- 
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observable" without further quahfication is one that is infrared safe at any scale, no matter 
how small. However, what is usually meant is that it is infrared safe at a scale not much 
smaller than the scale of the hardest interaction in the problem. What we want to do 
here is to specify the splitting scale at which the observable is sensitive to the splitting. 
That scale could be much smaller than Qq. 

Now we introduce the concept of the density operator p evaluated at resolution scale 
yu^. The idea is that interactions with scales greater than fs^ are included in |M)^M| while 
interactions with scales smaller than are integrated out (for final state interactions) or 
included in the parton distributions (for initial state interactions)."*^^ We can describe this 
in a rough way as follows. In each cut Feynman diagram, each integration region for final 
state partons that produces a collinear or soft divergence can be described as a region in 
which some group of partons with labels i become collinear to a given direction of a mother 
parton or some become soft, so that (X^jPi)^ 0. Divide this region into subregions with 
Cl2iPi)'^ < fj,'^ and (X^jPi)^ > In the (X^jPi)^ < region we can combine the partons i 
into a single effective parton for purposes of calculating the observable. That is, the partons 
i are "unresolved." This leads to a free integration over this region using a constant F for 
fewer partons. Adding these real emission integrals to the corresponding virtual diagrams 
and counterterms from the parton distributions gives a finite result containing logarithms 

of pVqI 

A more intuitive way of thinking about this is to imagine writing the Feynman diagrams 
in a coordinate space representation, in which we integrate over the positions xf of the 
interaction vertices relative to the position of the hard interaction, which is determined to 
within 1/Qg. Then we can restrict these integrations to \xf\ < 1/ p? before integrating over 
the final state momenta. 

In order that the parton distributions in Eq. (3.17) include initial state interactions at 
all scales smaller than fj?, the factorization scale //p at which the parton distributions are 
evaluated should be y?'. 

With this meaning of p evaluated at resolution scale /.i^, {F\p^ is invariant under 
p^ ^ p? + 5p'^ as long as F is infrared safe at a scale equal to p^ or larger. However, if 
one looks at with a resolution scale smaller than /i^ then (-^|p) does see the effect of 
changing p^ . 

Our object in this paper is to construct an approximate version of jp) as a function 
of the resolution scale. The idea is to construct jp) using a parton shower, starting from 
p^ = Qq and evolving to smaller values of /U^, down to a final infrared cutoff. This intuitively 
appealing idea has been inherent in the idea of a parton shower since the earliest days. 
It would be very useful to have a precise field theoretic definition of jp) as a function of 
resolution scale. However, this is beyond our scope in this paper. 

One finds that the logarithm of the resolution scale p'^ is more useful as a variable than 
p^ itself. Therefore we define Monte Carlo time t by 

/i^ = Qoe"* (3.37) 

^^This is the idea of the standard factorization theorem [20]. However, we here go beyond anything that 
has been proved. 
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and write the density operator with this resolution scale as \p{t)Y Then also the parton 
distributions are evaluated at factorization scale 

Ml = Qoe"* • (3.38) 

3.5 Parton shower evolution 

We are now prepared to set up a quite general framework for describing a parton shower. 
We take the framework as a set of axioms that we hope are reasonably intuitive. Later, we 
relate the operators that occur to the structure of Green functions at tree level in QCD. 

We use the evolution variable t that specifies the resolution scale of p according to 
Eq. (3.37). Thus t starts at zero and increasing t corresponds to decreasing virtuality. 
One stops evolution at an infrared cutoff t{ at which the use of an evolution based on 
perturbation theory is no longer appropriate. For instance, t{ might correspond to a 1 GeV^ 
virtuality. The evolving shower is represented by a state \p{t)) that begins with an initial 
state |/o(0)). The evolution is given by a linear operator U{t,t'), with 

\pit)) =Kit,t')\pit')) . (3.39) 

Here U{t,t) = 1. These operators have the group composition property 

U{t3,t2)U{t2,h)=U{t3,h) . (3.40) 

The class of evolution operators that will use is defined by two operators, 'Hi(t) and 
V(t), according to the differential equation 

j^U{t,t') = [ni{t)-V{t)]U{t,t') , (3.41) 

with initial condition U{t, t) = 1. 

The first operator, 'Hi(i), represents parton interactions and, in general, changes the 
number of partons and their momenta. We specify Tij^t) by giving its matrix elements 

{{pJ,s',c,s,c}m'\'Hi{t)\{p,f,s',c,s,c}m) ■ (3.42) 

In a lowest order shower, which we consider in this paper, the operator TCi{t) describes 
1^2 parton splitting, changing a state with m final state partons to one with m + 1 final 
state partons. 

The second operator that controls evolution, V{t), describes the effect of virtual graphs 
and the "unresolved" part of real emission graphs. In a lowest order shower, we do not 
account for the virtual graphs exactly, but rather account for only the infrared singular 
part of the virtual graphs, which can be deduced from the real emission graphs. That is, 
V{t) is determined from TCi{t). The operator V(t) does not change the number of partons 
or their flavors or spins, but can change their color states. 

We construct the shower algorithm in such a way that it conserves probability in a 
certain sense. If we were dealing with e~^ + e~ — > hadrons, we would demand that the devel- 
opment of the shower does not change the total cross section. For hadron-hadron collisions. 
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the total cross section does not have a weh defined perturbative expansion. However, if we 
start with a state \p{0)) that is "hard" in the sense that [{p, f, s', c', s, c}m|p(0)) is non-zero 
only for parton configurations with a large transverse energy, then we can demand that as 
this state evolves into the shower, the contribution from |/o(t)) to the total cross section 
does not change. The observable that measures the total cross section is 

Fi{{p,f}m) = l ■ (3.43) 

We will call the vector corresponding to -Fi simply (l|. Using Eq. (3.32), the inner product 
of (l| with a basis state is 

0-\{P,f,s',c',S,c}rr,) = {{s'}m.\{s}rn){{c}m\{c}rn) ■ (3.44) 

The statement that shower evolution leaves contributions to the total cross section 
invariant is 

{l\U{t',t)\p) = (l|p) (3.45) 

for any (suitably hard) state jp). This requirement leads to a relation between the matrix 
elements of V, and 7ii. To derive this relation we multiply Eq. (3.41) on the left by (l| 
and on the right by After using Eq. (3.45), we get 

= {l\[ni{t) -V{t)]U{t,t')\p) . (3.46) 

Since this holds for any suitably hard state |p) we have 

= (l|[?^i(t)-V(t)] . (3.47) 

We multiply on the right by \{p, f, s', c', s, c}m) to obtain 

{l\V{t)\{p,f,s',c',s,c}m) = {l\'Hi{t)\{p,f,s',c',s,c}m) • (3.48) 

At this point, we need to discuss the structure of the parton splitting operator TCi. In 
subsequent sections, we derive the form of TCi{t) based on the structure of QCD tree level 
matrix elements in the limit that two of the m+l partons become massless and collinear, one 
becomes massless and collinear with one of the beam directions, or one (a gluon) becomes 
soft. In this limit, the matrix elements take a factored form, hard {soft collinear). 
This factorization leads to the definition'^^ of 'Hi{t) in Sec. 9. At the moment, what we 
need is the structure of (l|7ii(i)|{p, /, s', c', s, c}m), which represents the inclusive splitting 
probability at splitting scale t. We will find 

{l\ni{t)\{p,f,s',c',S,c}m) =2{{s'yn^\{s}m){{c'}m\h{t,{p,fU)\{c}m) , (3.49) 

where the function h{t,{p, f}m) is given in Sec. 12, Eq. (12.21). The important point is 
that there is a trivial spin structure and a non-trivial color structure. 

^^The definition is not unique because there is freedom to choose what to do away from the soft and 
collinear limits 
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We take the operator V{t) to operate only on the color space and define its action in 
terms of its matrix elements, for which we use the notation ({c', c}m| V(t; {p, f}m)\W, c}m) • 
The definition of V{t) in terms of these matrix elements is 



V{t)\{p,f,s',c',S,c}m) = J2 |te/'S''C',S,c}m)({c',c}m|V(t;{p,/}m)|{c',c}„) 

{c',c}m 



(3.50) 



Thus, using Eq. (3.44), 



{l\V{t)\{p,f,s',c',s,cU)= {WU\{sU){{c'U\{c}m) 

{S',£}m (3.51) 
X {{c,c}m\V{t;{p,f}m)\{c,c}m) ■ 

If we insert Eq. (3.49) and Eq. (3.51) into Eq. (3.48) and cancel the spin factors, we get 

X] {{^'}rn\{c}ni) ({c , c}m\V{t; {p, f}m)\{c , c}„) 

{£',£}m (3.52) 
= 2{{c'}m\h{t,{p,fU)\{cU) . 

There is a simple way to satisfy this equation. We define 

{{c',cU\V{t,{pJU)\{c',cU) = J,{c}m.\h{t,{pJU)\{c}m) ^f']: 

+ {{^'}m\Kt,{pJU)\{d'U)^ ■ 



(3.53) 



What we have done here is to decompose {{c',c}m\y{t,{p, f}m)\{c',c}m) into two terms. 
In the first term, nothing happens on the bra side of the density matrix but there is a 
virtual correction on the ket side, while in the second term nothing happens on the ket 
side of the density matrix but there is a virtual correction on the bra side.^^ With this 
definition, when we perform the sums over {c}m and {c'}m in Eq. (3.52) using Eq. (3.23) 
we see that Eq. (3.52) is satisfied. 

We note that one could imagine solving Eq. (3.41) numerically in the form 

U{t + At, 0) |p(0)) = [1 - V{t)At] U{t, 0) |p(0)) + ni{t)AtU{t, 0) |p(0)) . (3.54) 

That is, one could use small time steps in which either one of the partons splits or else no 
parton splits and the weights for different color states are readjusted. However, this is not 
the way that shower evolution is typically constructed. 

To proceed down a more traditional path, we define two other operators, VE^t) and 
Vs{t) with sum 

VE{t) + Vs{t)=V{t) . (3.55) 



^®What V contains is the singular parts of the virtual corrections, which are related to the coUinear and 
soft singularities of the real emission diagrams, plus an "unresolved" contribution from the real emission 
diagrams. Thus we obtain V from Ti. 
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The distinction between VE(i) and Vs(t) lies in how we treat them within shower generation: 
VE(i) is exponentiated and Vs(t) is subtracted. We express the solution of Eq. (3.41) in 
the form ^ 

U{t,t')=M{t,t')+ [ dTU{t,T)[ni{T)-VsiT)]M{T,t') , (3.56) 

Jt' 

where the operator J\f{t,t') is the time ordered exponential of the operator VE(i), 

AA(t,t') =Texp|-^^drVE(T)| . (3.57) 

Here T represents the ordering in evolution time t. The operator M{t, t') is a generalization 
of the standard Sudakov exponential in parton shower Monte Carlo programs. It has the 
group multiplication property 

AA(t3, t2) Mit2,h) = Mit3,ti) (3.58) 
and satisfies the differential equation 

^^M{t,t') = -VE{t)Af{t,t') . (3.59) 

Eq. (3.56) is interpreted as saying that either the system evolves without splitting from t' 
to t, or else it evolves without splitting until an intermediate time r, splits or undergoes a 
color change at r, and then evolves (possibly with further splittings or color changes) from 
r to t. The first term contains a summation of effects from virtual splittings. In the second 
term we have a parton splitting contribution along with a subtraction that arises from the 
part of the virtual splitting contribution that was not summed to form part of Af{t,t'). 

Now the operator V{t) is completely defined but we still need to define VE{t) and Vs(i). 
Here we have some freedom. There are at least three obvious choices: 

1. We could define VE{t) = 0. Then Vs(t) = V(t). This choice leads to a trivial Sudakov 
exponential, J\f{t,t') = 1. This is similar to what one does in fixed order calculations 
when the singularities of the real emission graphs with m + 1 final state partons are 
removed by the subtraction terms with m partons. This is not useful in the context 
of a parton shower. 

2. We could define Vs{t) = 0. Then VE(i) = V(t). This means we exponentiate the 
whole virtual splitting operator. It is the most "shower way" to organize the parton 
evolution. The integral of V(t) over a range of t produces large logarithms and all 
of these logarithms appear in the Sudakov exponent. Since V{t) is a non-diagonal 
matrix in color space, the implementation of this choice may present difficulties. 

3. Alternatively we can define VE(t) to be the diagonal in color as follows: 

({C', C}^| VE(t, {P, f}m)\{c', CU) = C'fc 

X [{{cU\h{t,{p,fU)\{cU) (3.60) 

+ {{c'}m\h{t,{p,f}m)\{c}m)] • 
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Then Vs{t) is 



Vs{t, {p, fU) = V{t, {p, fU) - VE{t, {p, f}m) . (3.61) 

With this choice, the computation of the operator M{t, t') is simple because the basis 
vectors in the statistical space are eigenvectors of the operator V^it)- 

With this alternative, we do not exponentiate everything. However, when we study 
color in the following sections, we will see with this choice that (with the color basis 
that we will choose) Vs(t) is small compared to V-E.{t). First, we have 

{{c}m\h{t,{p,f}m)\{c}m)= ^ {{c}m\{c}ni) jl{c}rn\h{tApJ}rn)\{c}rn) , 

{c}m 

{{c}m\h{t, {p, f}m)\{c}m) = ^ ({c'}m | /l(t, {P, /}m) | {c}™)^ ({c}m|{c'}m) ■ 

{c}m 

(3.62) 

We will see in Eq. (7.14) that the matrix ({c}m|{c}m) is the unit matrix except for 
corrections. Thus the difference between using dual basis vectors |{c}m)^ in 
Eq. (3.53) and ordinary basis vectors |{c'}m) is not important in the large A'c limit. 
Second, we will see in Sec. 12 that the matrix {{c}m\h{t,{p, f}m)\{c}m) is almost 
diagonal in the sense that its off-diagonal matrix elements are suppressed compared 
to its diagonal matrix elements by factors of Thus the part of the virtual 

contribution that is not exponentiated is small. This small part is not neglected, but 
we can leave it out of the Sudakov exponent and treat it as a subtraction instead. 
Effectively, this means that we treat as a small parameter in addition to Og. 

4. Momentum and flavor mapping 

In this and the following sections, we explore how to define the splitting operator TCj. 
The first issue to examine is the momentum mapping. We begin with an m parton state 
with momenta {p}m- One of the partons, with label I S {a, b, 1, . . . , m}, splits. After the 
splitting, we have an m + 1 parton state with momenta {p}m+i- Our notation is that 
parton I splits into partons with labels / and m + 1, while the other partons keep their 
labels. The momenta {p}m+i after splitting are determined by the momenta {p}m and a 
momentum splitting variable that we call Cpi which defines the momenta of the daughter 
partons. There is also a flavor splitting variable, Q, which tells the daughter flavors. 

In this section, we first describe the splitting of a final state parton, then move on 
to the somewhat more complicated splitting of an initial state parton. For the final state 
splitting we first describe how {p,f}m+i is determined from {p,f}m and {CpiCf}- Then 
we state the inverse transformation, from {p, to {p, f}m and {CpjCf}- Finally, we 
deduce the jacobian for this change of variables. We will then be ready to do the same 
thing for the splitting of an initial state parton. 

There are many ways to define the momentum mapping {p}m ^ {p}m+i- One of the 
most successful is that of Catani and Seymour [14]. This may be called a local mapping: 
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the momenta of most of the partons are left unchanged, while the momenta of two partons 
are mapped into the momenta of three partons, {jpuPk) ^ ipuPm+iiPk)- Here parton k 
is a spectator parton, chosen as one of the partons that is color connected to parton /. 
(There is an exception to this rule in the case of initial state splittings.) The antenna 
factorization of Ref. [21] also uses a local mapping. We use a global mapping, in which all 
of the partons participate, as in Ref. [22]. This way, each parton has to contribute only a 
little momentum. We also include quark masses in the kinematics. 

4.1 Splitting a final state parton 

We begin by defining what happens to the parton flavors when a final state parton / splits. 
The partons with indices other than / and m + 1 keep their flavors, 

/i=/i , i^{/,m + l} . (4.1) 

What happens to partons / and m + 1 is given by the value of the variable = (//) fm+i)- 
The flavor splitting variable takes values in a set ^i{fi) that depends on the flavor of the 
mother parton. If parton I is a quark or antiquark, then the set ^i{fi) has only one element, 

^iUi) = {Uug)} , ///g • (4.2) 

Here, we have used the freedom to assign labels in order to assign the label I to the daughter 
quark or antiquark and the label m + 1 to the gluon. If parton / is a gluon, then <^f can be 
a pair of gluons or any choice of (g, q) flavors, 

<I>Kg) = {(g,g),(u,u),(d,d),...} . (4.3) 

In the case of a g ^ g + ^ splitting, we again use the freedom to assign labels in order to 
assign the label / to the daughter quark the label m + 1 to the daughter antiquark. 

We now turn to the momenta. Parton I has momentum pi with pf = m?'{fi) and splits 
into two partons / and m+1 with momenta and pm+i respectively. The daughter partons 
are on-sheh: pf = nT?{fi) and p^^j^^ = w?{fm+i)- We always have {pi + pm+i)"^ > m^i^fi). 

We need a bit of notation. Let Q be the total momentum of the final state partons, 

m 

Q = Y.P^ =Pa+Pb ■ (4.4) 

Define 



ai 



2prQ ' 



(4-5) 



2prQ 

Note that ai + bi > 1. To see this, let K = Q — pi be the total momentum of the final state 
spectator partons, with > 0. Then 

< 



2prQ 

-IprQ + m^fi) (4.6) 
2prQ 
ai + bi-l . 
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In order to define the momentum mapping, we first determine the total momentum 



Pl=Pl+Pm+l (4.7) 
of the daughters of parton /. We take Pi to be a linear combination of pi and Q, 

Pi = Xpi + ^"^^^ Q . (4.8) 
2a/ 

There are two parameters in this definition. The first, y, is a measure of the virtuality, 
P^ — m?{fi), of the splitting. The second. A, is a function of y that we will determine 
presently. 

For an exactly collinear splitting or the emission of gluon with momentum Pm+i = 0, 
we have Pi = pi- Away from these limits, the spectator partons will have to donate 
some momentum in order to allow Pi ^ pi. We elect to leave the momenta of the initial 
state partons unchanged, pa = Pa and ph = Ph- Instead, we choose to obtain the needed 
momentum from the final state spectator partons by letting the momenta after the splitting 
be related to the momenta before the splitting by a Lorentz transformation, 

p';=KPj , j^{l,m + l} . (4.9) 

With this method of transferring momentum, each parton donates a share of the needed 
momentum, with low momentum partons donating only a little momentum. 

The total momentum of the final state spectator partons before the splitting is 

K = Q-pi . (4.10) 

Since the momenta of the initial state partons remains the same, Q = Pa + Ph is the same 
as Q. The total momentum of the final state spectator partons after the splitting is then 

k = Q-Pi . (4.11) 

Since each final state spectator is changed by a Lorentz transformation, we have 

k^' = A'^^ . (4.12) 

In fact, there is a Lorentz transformation that does this, namely 

A{K, KY, -g>t ^K^Kf ' ^ ^ ^ 

provided that = K"^ . Thus Pi must lie on the hyperbola {Q — Pi)'^ = {Q — Pi)"^ in the 
Q-pi plane, as illustrated in Fig. 3. 

In the case in which the momenta K and K are carried by a single massless spectator, 
K is parallel to K. In this case, one can use an alternative representation of the boost in 
Eq. (4.13) that remains well defined when K'^ = {K + K)"^ = 0, 

A{k, KY, = < + - l) n^n, + n^n, , (4.14) 
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Figure 3: Momentum space mapping for a final state splitting. Coordinate axes in the two lightlike 
directions in the plane of pi and Q are shown. Points pi and Q are shown in the case that pi is 
nearly lightlike. Then P; lies on the hyperbola {Q — Pi)"^ = const, that passes through pi. 



where n and n are lightlike vectors in the Q-pi plane with n ■ n = 1 and {pi ■ n/pi ■ n) < 
{Q ■ n/Q ■ n). That is, these are the vectors along the two coordinate axes in Fig. 3, with 
n directed toward the upper right and h directed toward the upper left in the diagram. 
We can arrange that = by making a proper choice of A in Eq. (4.8). We have 

K'' = {ai + hi-l)2Q-pi , (4.15) 

and 

^2 = — [{2ai - 1 - y)2 - (1 - 4aibi)X^] 2Q-pi . (4.16) 
We obtain = if 



V 1 - 4ai6/ 

We note that A ^ 1 when y ^ 0. As y increases, A decreases. There is a maximum 
possible value of y, namely the value that makes A = 0, 

Wx = - Vai + k-iy -k . (4.18) 

The splitting parameter y has a simple meaning. From Eq. (4.8), we find, using 
Eq. (4.17), 

y = ' ^ ' ■ 4.19 

2pi-Q 

That is, y is a dimensionless measure of the virtuality of the splitting. An alternative 
relation is 

Pf-m\fi) 
2PrQ-{P^-m^{fi)) 
There are also alternative relations for ai and bi, 



2Pi.Q + w?Ui)-Pf ' 



(4.21) 

h 



2Pi-Q + m?{h)-Pf 
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These relations can be derived with the use of the relation 



2Pi-Q = {l + y)2pi-Q . (4.22) 

There is also a minimum value of y. Since Pf > [m{fi) + m{fm+i)]'^ , we have y > j/min 
where 

_ [m{fi) + m{fm+i)?-m''{fl) 

2prQ ■ ^ 

Note that if all of the partons are massless, then ymin = 0. In addition, if fi = // is a 
massive quark flavor and fm+i = then also ymin = 0. 

We have seen how, given the virtuality variable y, we can define Pi = pi + Pm+i in the 
plane of pi and Q such that P^ — m?{fi) = y2pi ■ Q and so that the needed momentum 
donation from the spectator partons can be obtained by a Lorentz transformation. It 
remains to define pi and Pm+i individually. This is simple. We let pi and Pm+i be any 
momenta on the appropriate mass shells that sum to P/. 

It will prove convenient to formulate this rather abstractly. We denote the daughter 
parton momenta by 

Cp = {Pl,Pm+i) ■ (4.24) 

We must have 

Cp^mP}m,Q) , (4.25) 
where the set /]({p}m) Cf) is described as follows, 

rii{p}m,Q) = [iPl,Pm+l) fi = m^ifi), PrQ > 0, 

Pm+l = '^^(/m+l), Pm+rQ>0, 

e,..af^{p'(+p';n+l)PlQ^ = 0, 

{Q - PI - Pm+lf = {Q - PI?} ■ 



(4.26) 



That is, each of pi and Pm+i lies on the appropriate forward mass shell, their sum lies in 
the plane of pi and Q, and = K'^. These conditions entail that 

[mCfl) + mUm+i)]'' < {pi + Pm+if < ( - \/Q'' + rn?{fi)-2pi-Q ) ' , (4.27) 

which corresponds to ymin < 2/ < 2/max- 

Wc note that the set Fi is a three dimensional surface in the space of momenta 
{PliPm+i)- One can choose three coordinates to describe this surface, for instance a virtu- 
ality variable, a momentum fraction variable, and an azimuthal angle. Different choices of 
coordinates may be best for different purposes, so we leave this choice open. 

We give a name to this transformation of momenta and fiavors: 

{p,/}„+i = i?K{p,/}m,{Cp,Cf}) • (4.28) 

In parton splittings close to the collincar or soft limits, partons lose energy when they 
split. It is of interest to see how this property carries over to splittings that are not close 
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to the limit. Here, the choice of the mapping Ri may be considered to be part of a model 
for shower evolution that could be sensible or perhaps not so sensible. An investigation of 
properties of the model is thus of some significance. 
Using = K^, we derive 

Q-pi-Q-pi=Pm+i-K+l[m\fl)-m\fi)+m\fm+i)] . (4.29) 

This is especially interesting in the case that parton I is a quark and parton m + 1 is a 
gluon. Then the term involving masses vanishes. We note that Pm+i ■ K > because 
both Pm+i and K lie inside or on the forward light cone. Thus the energy of the quark, 
as measured in the frame in which Q is at rest, is bigger before the splitting than after 
the splitting. By emitting bremsstrahlung, the quark slows down. If evolution of the final 
state were to continue long enough, the quark would slow to a stop in the Q = frame. 
Then radiation from that quark would cease. At any point in the shower evolution, there 
can be an initial state splitting, discussed below in Sec. 4.4. This changes pa. or and 
thus Pa +Ph = Q- Now final state quarks tend to come to rest in the new Q = frame. We 
judge that a tendency for quarks to slow down (and, similarly, for gluons to lose energy) 
is reasonably sensible. Of course, the partons should not be allowed to shower indefinitely. 
At some resolution scale, a perturbative model for showering is simply wrong and a process 
by which partons combine to form hadrons is needed. 

4.2 Combining two final state partons 

It is significant (and useful) that this transformation has an inverse. Let start with 
and determine {p}m and {Cpi Cf}- 

The splitting variable for the momenta is given by the momenta of the daughter par- 
tons, Cp = {pi,Pm+i)- From {p}m+i we determine 

■m+l 

'3 = 2^ Pi = Pa + Pb . (4.30) 

J=l 

Then Eq. (4.20) gives y, Eq. (4.21) gives ai and bi, and Eq. (4.17) gives A. Since the 
calculation of A involves taking the square root of A^, we should check that A^ > 0. For 
this purpose, we can express A^ in terms of dot products of vectors as 

4Q.k[Prk + P^-m^{fi)] + {P^-m^{fi)y ' 

Since Q, Pi, and K = Q — Pi lie inside or on the positive lightcone and P^ — m?'{fi) > 0, 
both the numerator and the denominator are non-negative. With y, A, ai and bi at hand, 
one can calculate the lightlike momentum pi by rearranging Eq. (4.8), 

PI = \{pi+Pm+i) - ^ " ^ Q . (4.32) 

We now have pi. We define p^ = Pa and pb = Ph- This leaves the pj for j ^ {/, a, b}. 
For this, we need the inverse Lorentz transformation to Eq. (4.9). From K = Q — pi and 
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K = Q — Pi — pm+i, we construct A{K,K)^i, using Eq. (4.13) or Eq. (4.14) with the roles 
of K and K interchanged. Then 

p>; = A{K, kr, , j i {I, a, b} . (4.33) 

The transformation of the flavors is simple. The splitting variable Cf is given by the 
flavors of the daughter particles, (//, /m+i)- The flavor of the mother parton is 

fi = fi + L+i , (4.34) 

with the obvious definition of adding flavors, as in d + g = d and u + u = g. The flavors of 
the other partons are unchanged 

fj = fj , j^{l,m + l} . (4.35) 

We give a name to this transformation of momenta and flavors, 

{fe/}m,{Cp,Cf}} = Q/({p,/Wi) • (4.36) 
This is the inverse transformation to Ri, Eq. (4.28). 
4.3 The integration measure for final state splitting 

With a suitable choice of the integration measure d(p for integrating over the splitting 
variables Cpi we can arrange that 



[d{p,f}m+i] g{{p,f}m+i) 

= I [d{p, fU] Yl [ '^^P ^(^P ^ ^KMm, Cf)) 9{{P, f}m+l] 

for an arbitrary function g{{p, f}m+i)- The definition that we need is 

PrQ 



(4.37) 



dCp = dy 9{y^i^ <y < ymax) A • 



vr 



2n6^{pf - m\h)) 2.K{frn+i - m\U,)) (4.38) 



(27r)4 ^-"^^ (27r)^ 

X {2nf 6 (^pi + pm+i - Xpi - ^ ^ Q 

Here the limits on y are given in Eqs. (4.18) and (4.23). 
4.4 Splitting an initial state parton 

Consider the splitting of an initial state parton, say parton "a." The initial state parton 
with momentum pa ~ VaPA splits to produce a new initial state parton with momentum 
Pa ^ VaPA and a new final state parton with label m + 1 and momentum Pm+i- We are 
using the usual backwards evolution here, so that the evolution going forward in time is 
Pa — > Pa + Pm+i- In this subsection, we describe how {p, f}m+i is determined from {p, /} 



m 
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and splitting variables {CpiCf}- The splitting of the other initial state parton is described 
by the same formulas with a <-> b. 

We begin by defining what happens to the parton flavors, recalling our notation that 
for the initial state partons "a" and "b", /a, fa, fh and /b denote the opposite of the flavor 
of the physical incoming parton. The partons with indices other than "a" and m + 1 keep 
their flavors, 

fj = fj , i^{a,m + l} . (4.39) 

What happens to partons "a" and m+1 is given by the value of the variable Cf = (fa, fm+i)- 
The flavor splitting variable takes values in a set <^a(/a) that depends on the flavor of the 
mother parton. This set is determined by the requirement that /a = /a + fm+i- If parton 
"a" is a quark or antiquark, then the set $a(/a) has only two elements, 

^a(/a) = {(/a,g),(g,/a)} , /a / g • (4.40) 

If parton "a" is a gluon, then Q can be a pair of gluons or any choice of {q, q) or (g, q) 
flavors, 

^a(g) = {(g,g),(u,u),(u,u),(d,d),(d,d),...} . (4.41) 

We now turn to the momenta. Let parton "a" with momentum fraction ?7a radiate a 
parton m + 1 with momentum Pm+i- We then need to define how to determine {p}m+i 
from {p}m and Pm+i- 

As discussed in Sec. 3.1, we can include masses for the initial state partons.^'' We take 
the partons to be on-shell with zero transverse momenta, 

m'^ifa) 

Pa = riaPA H P-B , 

VaS 

Ph = VhPB+'^^^PA , (4.42) 
m^fa) 

Pa = VaPA H : PB ■ 

VaS 

Recall that we define pA and pB to be lightlike approximations to the incoming hadron 
momenta, with 2pA ' PB = s. The radiated parton can have a mass, 

Pl+i = m\fm+i) ■ (4.43) 

We take the momentum fraction of parton "b" to remain the same, 

flh = Vh ■ (4.44) 

The momentum fraction fj.^ after the splitting will be determined by Pm+i- As in the case 
of a final state splitting, it is not generally possible to have Pa = Pa + Pm+i given the 
mass shell conditions and the possibility that the radiated parton has non-zero transverse 



^^Recall that this is optional. One could just replace the masses for the initial state partons by zero, 
^^(/a) = TTT-^ifa.) = {fh) ~ 0. One could also set all quark masses to zero. 
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momentum. In order to allow the approximation that both p^. and pa are on-shell with 
zero transverse momenta, we therefore take some momenta from the final state spectator 
partons by letting the momenta after the splitting be related to the momenta before the 
splitting by a Lorentz transformation, 

p^^=K^'^p- , j G {!,..., m} . (4.45) 

Since each final state spectator is changed by a Lorentz transformation, we have 

k^' = A'^^ K"" , (4.46) 

where K is the momentum of the final state partons before the splitting, 

K = p^ + p^ , (4.47) 

and K is the momentum of the final state spectators after the splitting, 

K = Pi, + Ph - Pm+i ■ (4.48) 

In order for K and K to be related by a Lorentz transformation, we need = K^. 
To see what this means, define 

Q(rya) = Pa + Pb = (f,. + PA + f % + PB ■ (4.49) 

Then we demand that 

{Q{f,^) - Prn+l? = K"" . (4.50) 

The vector Pm+i determines r)a. In the space of Pm+i^ a surface of constant ?7a is the 
intersection of the hyperbola = rn'^ifm+i) with the hyperbola given by Eq. (4.50), as 

illustrated in Fig. 4. Using = rnPifm+i) hi Eq. (4.50), we can write 

2Q{fj,)-p^+i = Q{f,,f + m\fm+i)-K^ ■ (4.51) 

Looked at this way, a surface of constant r/a is the intersection of the hyperbola p^+i — 
m^(/m+i) with the plane defined by Eq. (4.51). We will require rja < 1. Thus the allowed 
region in pf^^i is the part of the forward mass shell with 

2Q{l)-Prn+i<Q{lf + 7n\fm+i)-K^ . (4.52) 

In order to solve for r/a given Pm+i-, we write and in the form 

2 P 

K = aijs, 7 , 

K = afji,- 7 . 



(4.53) 
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Here 



while 



P=- ^"^f^)^"^f^) , (4.54) 



a = r]i,s -2pA-Pm+i , 

P= i2p-B-Pm+i } , (4-55 j 



The condition K = K now determines r/a, 

1 



r/a = ^|i^'+7+ V(^'+^)'+4«/3| • (4.56) 
It is a consequence of Eq. (4.56) and the kinematic conditions (3.6) and (3.9) that 

??a > % • (4.57) 

We prove this in Appendix A. 

Having fixed r)a so that = K^, these two momenta will be related by Eq. (4.46), 

k^" = A{k, K)% K" , (4.58) 

where A^i, the Lorentz transformation (4.13). This allows us to define the spectator mo- 
menta after the splitting to be related to the spectator momenta before the splitting by 
this same boost, as in Eq. (4.45), 

pi; = A{k, KY.p'j , i G {1, . . . , m} . (4.59) 

For a final state splitting, we defined a splitting variable Cp in a three dimensional set 
Fi so that {p}m together with determines {p}m+i- Here we use an analogous notation, 
with 

Cp = (Pa,Pm+l) • (4.60) 

We choose 

Cpe/^a(Mm,Cf) , (4.61) 

where 

A(Mm,Cf) = |(^a,Pm+l) pli+l = "^^(/m+l), An+r(Pa +Pb) > 0, 

Pa = r?aPA + — PB, 

??aS (4.62) 
k^{n,)=K^{r]^), k-K>0, 

r?a < -^a < l| . 
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Figure 4: Momentum space for Pm+i in an initial state splitting. Coordinate axes along pA and 
Pb are shown, along with the point Q(?7a). We show the mass shell for Pm+i and the hyperbola 
(QiVa.) ^ Pm+i)^ = K^. Points Pm+i for a given value of 77a are on the intersection of these two 
hyperbolas. We illustrate a point Pm+i with zero momentum transverse to p^ and pb, but in 
general Pm+i can have transverse components. The plane specified in Eq. (4.51) is also depicted. 

That is, the radiated parton lies on the forward mass-shell and the new initial state parton 
is also on its mass shell with has zero transverse momentum. Its momentum fraction f)a is 
determined by = K^, with K inside the forward lightcone. The solution of = 
is given by Eq. (4.56). We also need 7)a to be big enough that that the kinematic bound 
that was given in Eq. (3.9) is obeyed. This follows from 57a > r/a- Finally, r/a must be less 
than 1. 

As for a final state splitting, the set Fg. is a three dimensional surface in the space 
of momenta (pa,Pm+i)- One can describe this surface using three coordinates such as a 
virtuality variable, a momentum fraction variable, and an azimuthal angle. We leave the 
choice of coordinates open. 

We give a name to this transformation of momenta and flavors, Ri with / = a: 

{p,/Wl = i?a(fe/}m,{Cp,Cf}) • (4.63) 

4.5 Combining an initial state parton with a final state parton 

The transformation for splitting an initial state parton has an inverse. Let us start with 
{P,f]m+i and determine {p,f]m and {Cp,Cf}- 

The momentum splitting variable is simply i^p = (pa,i^m+i)- For the momentum frac- 
tion of parton "b" before the splitting, we take 

r?b = r)b • (4.64) 

To determine r/a, we simply use the representation (4.53) of and and solve = 
for r^a in terms of fja.. This gives 

r/a = z^|i^'+7+ V(^'+7)'+4a/3| . (4.65) 
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Once we have r^a, we can construct K = Pa,+Ph- Then from K and K we can construct 
A{K,K), the inverse Lorentz transformation matrix to A(K,K) defined in Eq. (4.13) by 
simply using the same equation with K ^ K. Then we can construct the remaining 
momenta using 

p>; = A{K,kr,p'^^ . (4.66) 

The transformation of the flavors is simple. The splitting variable (f is given by the 
flavors of the (backwards evolution) daughter particles, (/a, fm+i)- The flavor of the mother 
parton is 

h = h + L+i ■ (4.67) 
The flavors of the other partons are unchanged 

, j^{a,m + l} . (4.68) 

We give a name to this transformation of momenta and flavors, 

{{pJU,{Cp,Ci}} = Q.{{pJ}m+i) ■ (4.69) 

This is the inverse transformation to i?a, Eq. (4.63). 

4.6 The integration measure for initial state splitting 

With a suitable choice of the integration measure d^p for integrating over the splitting 
variables <^p, we can arrange that 

f}m+l] g{{P, f}m+l) 

= / [d{p, f}m] Yl / '^Cp ^(Cp e r,{{pU, Cf)) 9{{p, fU+i) 

for an arbitrary function g{{p, f}m+i) with support in > 4m^. The deflnition that we 
need is 

dCp = ^TT^ 27r5+(p„+i - m (/„+i)) ^— . (4.71) 

The factor {a+ [5 / rj^) / {a+ 13 / fj^) here is just dfia/drji^^ calculated from the relation = K^. 
5. Spin states 

The quantum scattering amplitude |Af({p, /}„)) is a vector in spin® color space. Thus 
we can expand it in terms of spin and color basis vectors, 

|M({p, /},„)> = |{c}^> Y K^}-> ^^^{P^ /' ^' c}-) • (5-1) 

The treatment of color that appears to us to be most useful for parton showers is a bit subtle 
and, in particular, involves basis vectors that are not exactly conventionally normalized 
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and are not exactly orthogonal to one another. The color basis is described in Sec. 7. In 
contrast, our spin basis vectors are quite standard, are orthogonal and normalized. 



({s'}m|{s}m) = <^m',m<^{s'}„,{s},„ • (5-2) 

The spin labels {sa; Sb, si, . . . , Sm} represent the helicities of the corresponding particles. 

The basis states for parton I are represented for quarks or antiquarks by Dirac spinors 
U{pi,si) or V{pi,si), where = m?{fi). For gluons, we need polarization vectors e{pi,si) 
with pf = 0. Our definition makes use an auxiliary vector ni, chosen along the intersection 
of the positive lightcone with the plane of pi and Q, where Q = YlY=i P'J ~ Pa + Pb is the 
total momentum of the final state particles. The normalization of n; is not important. We 
take the solution that is not close to pi in the case that m'^{fi) is small. A convenient set 
of choices is 



PB , 1 = a , 

PA , l = h , 

Q = pi , /G {!,..., ml 

Q-Pi + ^/{Q-PiY-Q^m^fi) 



For quarks and antiquarks, we use Dirac spinors U{p, s) or V{p, s) with p^ = m? and 
s = ±1/2. We can take V{p, s) = [U{p, s)C]^, where C is the charge conjugation matrix, 
27^7" in the chiral representation of the gamma matrices, with C~^^^C = —{'J^)'^ and 
= = C'^ = —C. The spinors obey {f — m)U{p,s) = and {p + 'm)V{p,s) = 0. 
They are normalized to 



[/(p,s)7^^7(p,s) =2p^ , 
V{p,s)j^'V{p,s)=2p'' . 



(5.4) 



We use helicity eigenstates, defined so that 



75^t/(p, ±1/2) = ±[/(p, ±1/2) , 
j5^V{p,±l/2) = ±V{p,±l/2) . 



(5.5) 



Here the spin vector s is 

1 m 

s = — p n , (5-6) 

m p-n 

where n is the auxiliary lightlike vector from Eq. (5.3). Thus = —1 and s • p = 0. 

A convenient definition that defines the phase of U{p, —s) in terms of the phase of 
U {p, s) is 

m 



U{p,-s) = 1± i^]V{p,s) . (5.7) 



^*Many authors follow the convention that the helicity label for an incoming particle is the negative of 
the particle's helicity. In contrast, our convention is that Sa and Sb, as well as si, . . . , Sm, represent the 
physical helicities of the particles. 
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For gluons, we need polarization vectors £^{p, s; Q), representing a given helicity s and 
defined with the aid of an auxiliary vector Q. The polarization vectors obey p-e{p, s;Q) = 
and 

Q-e{p,s;Q) = . (5.8) 

We can also write 

e{p,s;Q) = e{p,s;n) , (5.9) 

where n ■ e{p, s;n) = and n is a lightlike vector defined in Eq. (5.3). We define the phase 
by using the standard definition in terms of mass zero Dirac spinors [23] , 

.>'(p.±l;.) = ±^''!!^^/^'^''^"'-^'/^' . (5A0) 
\/2C/(p.Tl/2)C/{n,±l/2) 

With this definition, there is a simple relation between polarization vectors defined 
with different auxiliary vectors Q and thus different lightlike auxiliary vectors n [23], 

e^ip ±1- n) - e^ip ±1- n') = T= ^2 17(n, ±l/2Mn-, ^1/2) 

^ ' ^ ^C/(n, ±l/2)C/(p, Tl/2) C/(p,±l/2)[/(n',Tl/2) 

(5.11) 

If we use these polarization vectors with the exact tree-level Feynman diagrams, gauge 
invariance of the matrix elements, together with the fact that £^{p, s; n) differs from 
e^{p,s;n') by a vector proportional to p^, shows that the amplitude is independent of 
the choice of the auxiliary vector Q. Our matrix elements will be approximate and will be 
gauge invariant only to the extent that the splittings are close to the soft or collinear limits. 
Thus some dependence on the auxiliary vector used to define the polarization vectors will 
result. 



6. Splitting functions for the quantum states 

Consider for a moment a theory without spin, color, or flavors, say (p^ theory in six di- 
mensions as in Sec. 2. An (m -|- l)-parton scattering amplitude M({p}fn+i) is simple in 
the limit in which two of the partons are approximately collinear. Supposing that partons 
m + 1 and / are almost collinear, we have 

M{{pU+i)^vi{pU+i)M{{p}m) , (6.1) 

where v{{p}m+i) = g/C^Pi • Pm+i)- Here {p}m is determined from {p}m+i by the (six- 
dimensional version of) the transformation described in Sec. 4. This factorization formula 
for the amplitude becomes exact in the limit that partons m + 1 and I become collinear. 
Away from the collinear limit, there is some freedom to choose the momentum mapping 
and the splitting amplitude v{{p}m+i)- One has to make a definite choice based on ease of 
computation or conceptual simplicity. In the case of QCD, we have soft as well as collinear 
singularities, we have parton flavors (which are rather trivially treated) and we have color 
and spin, which are not so trivial. Let us see how to describe splitting in QCD. 
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It has been known for a long time that QCD amplitudes factor in the soft and coUinear 
limits [24]. Indeed, there are beautiful modern formulas for the factors [23, 25, 21]. We 
have adopted a more pedestrian approach that has at least the advantage of encompassing 
the soft and collinear limits at the same time and of including masses. A treatment of the 
squared amplitude that is rather similar to the approach of this paper, but at higher order, 
may be found in Ref. [26]. 

6.1 Definition of the splitting functions vi 

The QCD scattering amplitude for m + 1 partons is a vector |M({p, /}m,+i)) in color (^spin 
space. In the limit that two partons, / and m + 1 are almost collinear, this amplitude takes 
a certain limiting form, 

\M{{p, /}™+i)> ~ \Mi{{p, , (6.2) 

where |M;({p, f}m+i)) is to be defined precisely below. When Pm+i becomes soft, then all 
of the |M;({]5, /}m+i)) amplitudes contribute to the limit, 

|M({p,/Wi)>~El^^'({^'/Wi)> • (6.3) 
I 

We arrange the definition so that Eqs. (6.2) and (6.3) are exact in the collinear or soft limit 
respectively. We also arrange that |M/({p, /}m+i)) is defined for any {p, f}m+i- Then these 
equations are approximate away from the limit. The amplitude | Af;({p, /j^+i)) is then 
the contribution to the (m + l)-parton amplitude from the splitting of parton / in the 
parton shower approximation. 

We now need to define \Mi{{p, f}m.+i))- This amplitude factors into a splitting op- 
erator times the m-parton matrix element evaluated at momenta and flavors {p, f}m de- 
termined from {p,f}m.+i according to the transformation Qi{{p, f}m+i), Eq- (4.36) or 
Eq. (4.69), 

\Mi{{p, fU+l)) = tjifi ^fi + fm+l) ^/({P, /}™+l) \M{{p, f}m)) . (6.4) 

In Eq. (6.4), V"^ {{p, f}m+i) is the analogue of v{{p}m+i) but is now an operator on the 
spin part of the color® spin space. There is also an operator tj{fi fi + fm.+i) on the 
color part of the color fS> spin space. This operator multiplies by the right color matrix. 
We will not comment further on it in this section, but will turn to the description of color 
in Sec. 7. The spin dependent splitting operator can be described in terms of its matrix 
elements, 

\Vi\{p, f}ra+l)\{s}rn) ■ (6.5) 

This is a simple function of {p, f}m+i, {s}m+i, and {s}m- Furthermore, we can take it to 
be diagonal in the spectator spins, 

m+l = JJ 5s^^sA Vl{{p,f}m+l,Sm+l,Sl,Sl) . (6.6) 

\j<^{l,m+l} J 

Our object in this section is to define the splitting functions v from the QCD vertices. In 
writing formulas for vi{{p, f}m+i, Sm+i, si, si), we will use the momentum pi. We under- 
stand that this is obtained from {p, f}m+i according to the transformation Qi{{p, f}m+i)- 
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6.2 Initial state q ^ q + g splitting, quark scatters 




Figure 5: Illustration of factorization at the amplitude level, leading to the definition of the 
splitting function in Eq. (6.19). When the final state gluon becomes coUinear with the initial state 
quark, the full amplitude is approximately the amplitude for one fewer parton, with the incoming 
quark after the gluon radiation approximated as being on shell, times a splitting function simply 
made from the QCD qqg vertex and tlie singular quark propagator. If the final state gluon is soft, 
then this diagram is one of the possibilities. Then, the soft gluon could have been emitted from 
any of the external legs of the diagram and we must sum over all of the possible emissions. In this 
case, we can, however, use a simpler splitting function. 
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= e/i(Pm+l; Sm+lj Q) 



Consider an initial state g — > g + g splitting in which the gluon goes into the final state 
and the quark enters the hard scattering, as illustrated in Fig. 5. We highlight this case 
because it exhibits some complications compared to final state g' — > g + g splitting. The 
parton that splits could be either of the initial state partons. We examine the case that 
it is parton "a." The kinematics were defined in Sees. 4.4 and 4.5. In order to keep our 
notation for this subsection as simple as possible, we will write 



(6.7) 



Here e is defined to be orthogonal to Q = Pa + Pb = Pa + Pb- 

The Feynman rules for |M({p, /}m,+i)) give the following structure 

P A- m 

^ = ^ L , gtr^(p,s) . (6.8) 

Here there are factors ?7(p, s) for the initial state quark and e* for the final state gluon. 
There is a vertex (77^ with a color matrix . There is a propagator for the off-shell quark 
that carries momentum 

P=p-q . (6.9) 

The quark propagator has a denominator — m? and a numerator f + m. Finally, there 
is the rest of the diagram, H, which contains hard interactions. Thus H carries a Dirac 
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spinor index. We need to relate this to \M{{p, f}m)), where the momenta and flavors 
{Pi /}m as well as the splitting variables {CpiCf} are given by the momentum and flavor 
mapping Qi{{p, f}m+i), Eq. (4.36). 

We will be concerned with the behavior of M in the soft and coUinear limits. The soft 
limit is (7 — > (for all four components of q). Since — w? = — 2p • g, the denominator is 
proportional to a single power of g, while the numerator stays finite in the g — > limit. To 
describe the coUinear limit, let qx_ be the part of q orthogonal to pA and -p-Q. The coUinear 
limit is g_L — > while q • pB stays finite. Also, m ^ with rr? < Then near the 

coUinear limit, the denominator is 

P''-m'={J^ql-^-^m^{\^0{ql,m')) . (6.10) 
\ri-ri 7] J 

That is, the denominator has two powers of q±. To analyze the numerator, we note that 

p = 'qpx + 0{q\,m?) , 

p = fjpA + 0{ql,rr?) , (6.11) 

q = {'>)-Ti)PA + qL + 0{q\,w?) . 

The numerator has a factor [f + m)^*U{p, s). Using q - e = and {p — m))U{p, s) = 0, one 
can rewrite this factor in the form 



+ m)fU{p, s) = 2(p- -J^q) ■ e* U{p, s) + 
\ rj — r] / 



[jp — ml 

7] 



Uip,s) . (6.12) 



Looking at the numerator in this form, we see that it vanishes in the coUinear limit pro- 
portionally to one power of q± or m. Thus M is only half as singular in the coUinear limit 
as it first appears. 

We now develop an approximation for M. We can insert a factor 1 next to H in 
Eq. (6.8) so that it reads 



P2_^2 ^9t i )U{p, s) . (6.13) 

Here n is the lightlike vector n = pB • 

We now notice that the contribution from the first term, namely 



M,, = H^ -^2^ {gtr)U{p,s) , (6.14) 

can be neglected because it does not have a coUinear or soft singularity. To see this takes 
a little analysis. First, we write 

p = P+{p + q-p) ■ (6.15) 
Since {fi — m){f + m) = P^ — m^, we have 



M^, = H-^ {gty*)U{p,s) + H {gtr)Uip,s) . (6.16) 

2n-p 2p-n P'^ — 
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The first term is non-singular because the denominator is cancelled. In the second term, the 
vector {p + q — p) vanishes in the collinear or soft limit. In the soft limit, it is proportional 
to one power of q, which cancels the single power of q in the denominator. In the collinear 
limit, it is proportional to one power of q±. As we have just seen, the rest of the numerator 
contains an additional factor of q±. Together, these cancel the two powers of q± from the 
denominator. Thus no singularity remains. 
We are left with 

+ m)ii P + m 

Msing = , {9tr)U{p, s) . (6.17) 

In Eq. (6.17), the factor p + m is 

f + m = ^U{p,s)U(p,s) . (6.18) 

s 

The factor f/(p, s) is to be associated with H, giving the hard scattering amplitude for 
an incoming quark with spin s. The remaining factor, U{p,s), then becomes part of the 
splitting function. This calculation leads us to define the splitting function as 

V4vras ^ U{p,s)'^{p- ^ + m)'y^'U{p,s) 



(p — q)"^ — vn? ^ 2p-n 

This does not include the color matrix and a factor —1, which will be included in the color 
operator tj{fi fi + fm+i)- As indicated by the derivation, Eq. (6.19) is directly given by 
the factorized structure of QCD Feynman graphs in the soft and collinear limits. There is 
freedom to choose the form of the splitting function as one moves away from these limits. 
We have made a simple choice. 

In the hard part of the diagram, we can make approximations that are valid for — > 0. 
In particular, we can adjust the momenta of the partons with indices other than "a," 
replacing pj by pj as defined by the momentum mapping Ra.{{p, f}m, {Cp; Cf})) Eq. (4.63). 

6.3 Initial state q ^ q + g splitting, gluon scatters 

We consider next the process in which an initial state quark with label "a" splits to make 
a quark that goes into the final state (with label m + 1) and a gluon that enters the hard 
scattering (with label "a"). We simplify the notation as in the previous subsection by using 



(6.20) 



The Feynman rules for |M({p, /}m+i)) give the following structure 

n ( p- n) 

M = H^ U{q, s'){9t'Y)U{p, s) . (6.21) 
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Here H is the hard part of the graph, now with a vector index, and there is a propagator 
for the off-shell gluon that carries momentum 



P=p-q . (6.22) 

We have chosen the axial gauge n ■ A = 0, where n is the lightlike vector n = Pb- The 
numerator of the gluon propagator is 

D^^P; n) = -g'^-' + ^ . (6.23) 

As in the previous subsection, M is singular in the collinear limit, in which q± ^ 
and m — > 0.^^ As in the previous section, the numerator is proportional to q± and the 
denominator is proportional to q'^ in the collinear limit. 

In order to find a suitable approximation for M, we insert 1 next to H, so that it reads 



M = Hn 



-D°''^{p;n) + 



p-n 

We can drop the term {p°^n^ + n°'p^) / {n-p). To see this, we write 
p'^ni^ + n'^pf' P^n" + n"'P^ 



n (p- n) — 

p2 U{q,s'){gt'r)U{p,s) . (6.24) 



p-n P-n 



+ D'"'{p;n)-D''f'{P;n) . (6.25) 



The term P'^n^ gives zero when contracted with -D^u^. When we contract n°'P'^ with Dfj,,^, 
we get 

n"P^ Di,u{P\n) ^ n'^n, 

P-n P2 (p.^)2 ' i^-^t'i 

which does not have a collinear singularity. Finally, the difference D°'^(p;n) — W^^^^P-^n) 
is proportional to q_\_ in the collinear limit because P — p oc g_L in this limit. The remaining 
numerator factor gives another factor q±^ in the collinear limit, so that the factor q\ from 
the denominator is cancelled. 

Thus we are left with M ~ Mging, where 

M,i,g = -H^ D'-^ip; n) ^^^-^^^ U{q, s'){gt'inU{p, s) . (6.27) 
The factor D'^^{p; n) is 

D"'^(p;n) = 5;e"(p,s;n)£^(p,s;n)* = ^£-(p,5;Q)£^(p,s;g)* . (6.28) 

s s 

It is equivalent to use n = or Q = Pa + Pb = Pa + Pb to define the polarization vectors 
since they are orthogonal to both p = Pa and pb. The factor £"(p, s;Q) is to be associated 
with H, giving the hard scattering amplitude for a gluon with spin s. The remaining factor, 
e^{p, s; Q)*, then becomes part of the splitting function. 



^®We do not need to be concerned with the soft hmit, g — > with m —^ 0. Here, there is only a 1/^^ 
singularity, which is too weak to create a logarithmically divergent integration over final states. Nevertheless, 
our approximation to M also matches the behavior of M in the soft limit. 
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This calculation leads us to define the splitting function as 



"^' ■e^^ip, s; QT D,,{p - q; n) U{q, s')YU{p, s) . (6.29) 



iP - qf 



This does not include the color matrix, which is included in the color operator t\{fi 
fi + /m+l)- 

In the hard part of the diagram, we can again make approximations that are valid 
for q'j_ — > 0. In particular, we can adjust the momenta of the partons with indices other 
than "a," replacing pj by pj as defined by the momentum mapping -Ra({p, /}m, {Cp; Cf}); 
Eq. (4.63). We will also want to use Q = Pa + Ph instead of Q = Pa + Pb to define the 
polarization vector for parton / as it enters the hard scattering. This is just a change in 
notation rather than an approximation since Q is in the p-Q plane. 



6.4 Other qqg splittings 

For splittings involving a qqg vertex, there are several other choices for the flavors //, /; 
and fm+i and for whether the index / of the parton that splits is a final state index or 
an initial state index. The results for vi{{p, f}m+i, Sm+i, si, si) are listed in Table 1. In 
constructing this table, we keep track of two sign factors. First, there is a numerator sign 
that we compute as in the examples in the preceding subsections. This is a minus sign 
whenever a gluon or an antiquark enters the hard scattering from the initial state or leaves 
the hard scattering to the final state and a plus sign when a quark enters or leaves the hard 
scattering. Second, there is a color sign. There is always a color matrix f^. We count a color 
3 line entering the hard scattering as a color 3 line leaving the hard scattering. Viewed this 
way, the color matrix is actually {t"')'^- However, the generator of color rotations for the 
3 representation of SU(3) is —{t"')^. Thus we include — (t°)^ as part of the color operator 
tlifi —^fi + fm+i) and include the minus sign as part of the splitting function vi. This 
leaves a factor —1 in vi whenever a 3 line leaves the hard scattering. The sign included in 
t]{fi ~^ fi ~\~ fm+i) is indicated in the last column of Table 1. The net sign of vi is shown 
in the table. 

The construction for a final state splitting, as for an initial state splitting, makes use 
of a lightlike vector n/ that is in the plane of pi and Q = Pa + Pbi which is the same as the 
plane of and Q = Ps,+Ph- The normalization of n/ is not significant. Our choice for this 
vector was defined in Eq. (5.3). 



6.5 Splitting with a ggg vertex 

We construct the splitting function for a g — > g + g splitting in a similar fashion. In the 
case of a final state splitting, we use the ggg QCD vertex, 

v'^^^iPa^Pb.Pc) = g'^^iPa - PbV + /^(P6 - PcT + 5^"(Pc " Paf , (6.30) 
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Q) 



, U{pu k)l^\Pi + ^^+1 + mUi)]iiU{pu si) 



2pi-ni [{pi+Pm+iY - m?{fl)] 

^ yipi,si)ii[Pi - ^ifi)h^yipi,si) 

2pi-ni [{pi +pm+iY - m^ifi)] 

* U{PU sQjiiiPi - p^+i + m{fi))-f''U{pi, si) 
2pi-ni [{pi -pra+iY -m?{fi)] 



. A^* ^^i)^^^i>i - i'ni+i - m{fi))iiiV{pu si) 

-£fi{Pm+l,Sm+l',Q) — 



q q 



-e^{puSi]Q)D^y{pi +pm+i,ni 
-e^\pi, si; Q)*D^u{pi - pm+i; ni) 



2pi-ni [{pi -pm+iY -m2(/;)] 

U {pu Sl)YV{pm+l,Sm+l) 



{Pl +Pm+lY 

U(pm+l,Sm+l)YU {PU Si) 
{pi -Pm+i)2 



q q 



uf A\*r> f - ^,y{PhSl)Yy{Pm+l,STn+l) 
-£^{Pl,si;Q) D^y{pi-pm+l]ni) , 

\Pl - Pm+l) 



-£:,i{pi,si;Q) 



-^iiiPhSi^Q) 



U{pm+l,Sm+l)l^\Pl -f>. 



m+l 



m{fi)]'ifiV{pi,Si) 



2pi-ni [{pi -pm+i)'^ - frfi(fi)] 
U{pusi)rli\pi-p^j^^ + m{fi)\-i^'V{pm+i,Sm+i] 



2pi-ni [{pi 



m^Ul)] 



Table 1: Splitting functions vi{{p, /}m+i, Sm+i, si, si) for splittings involving a qqg vertex, with a 
common factor ^/Anag removed. The values of / are either in the set of initial state indices / = {a, b} 
or in the set of final state indices F = {1, . . . , m}. The flavors and fm+i can be / = g or can 

be a quark index, q € {u, d, . . . }, or an antiquark index q g {u, d, . . . }. Recall from Eq. (3.2) that 
in the case of an initial state parton, /( and // denote the opposite of the incoming flavors of the 
parton. For a row in which q denotes a quark index, q denotes the corresponding antiquark index. 
The next column gives the value of vi corresponding to the values of Z, and fm+i indicated. 

The last column indicates the sign of the color matrix that is incorporated into the color operator 
tj. The lightlike vector ni is defined in Eq. (5.3). 



to define 

Vl{{p, f}m.+l,Sm+l, Sl, Sl 



, ea(Pm+i, Wi;Q)*e/3(Pi>s«;0)*e''(P/,s«;Q) (6.31) 
2pm+i -pi 

X V°-^^ {pm+\,Ph -Pm+1 -Pl) D^uiPl +Pm+i;ni) . 



- 44 - 



For an initial state splitting, we have 

Vl{{p, f}m+l,Sm+l,Sl,Sl) 

= - ea{Pm+l, Sm+i; QT ^ 01, k, {Ph Si; QT (6.32) 

-Pl 

X V°''^^{pm+l,-pi,Pl - Pm+l) D^viPl -Pm+i;ni) . 

In each case, we have the exact QCD vertex and the exact propagator for the ofF-shell gluon 
in • ^ = gauge followed by a projection onto the physical gluon degrees of freedom 
contained in the on-shell polarization vector. 

6.6 Soft splitting function 

These splitting functions enable us to approximate the (m + l)-parton matrix element in 
the cases that Pm+i is collinear with pi or else Pm+i is soft. In the special case that Pm+i is 
soft, or possibly soft and collinear with pi, a simpler splitting function can be used. When 
Pm+i is soft, we have 

\M{{p, fU+i)) ~ Yl \Mr'\{P, fU+i)) , (6.33) 
I 

where 

|Mrft({^, /}^+i)> = t]{fi -^fi + L+i) Vi^'^''{{p, f}m+i) \M{{p, fU)) . (6.34) 
The matrix elements of Vi'^°^^ are specified by a function vf"^, 

({S}^+l| V;^'^°"({^, f}m+l)\{s}m) = ( n I ^"''({P, f}m+l, Sm+1, k Si) . 

(6.35) 

If parton m + 1 is a quark or antiquark, f^"^* = 0. When parton m + 1 is a gluon, 

soft/r- -fl ^ ~ N. /-. r ^{Pm+ljSm+l'iQ)* 'Pl (a oa\ 

Vi {{p, /}m+i, Sm+1, su Si) = V47ras Ss,,si 7 7 ■ (6.36) 

Pm+l-Pl 

The functions are not as powerful as the functions vi because they provide good 
approximations only in the soft gluon limit. Nevertheless, we will have occasion to make 
use of them. 

7. Description of color 

We will need a description of the quantum color state that is adapted to a description 
of shower evolution. If we use an index notation, to each parton with label I there is 
associated a color index a/, which takes values 1, . . . , 3 for a quark or antiquark and takes 
values 1, . . . , 8 for a gluon. There is also a spin index A^, which takes values ±1 /2 for quark 
and ±1 for a gluon. We can expand M in terms of color and spin basis vectors in the form 

{cU {S}m (7.1) 

X M{{pJ,S,c}m) , 
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where the ^{{c}m) form a basis for the space of color singlet amplitudes with color labels 
{c}m and the E{{s}m) form a basis for the spin space with spin labels {s}m- The quantities 
M{{p, /, s, c}) are the expansion coefficients. In a vector notation, this is 

Here |-/Vf({p, f}m))^ g lies in the combined color-spin space while |{c}m)j, is a vector in color 
space and |{s}m)g is a vector in spin space. 

As discussed in Sec. 5, we use a conventional treatment of spin. We assume that the 
spin labels A already represent parton helicities using suitable conventions for choosing 
corresponding Dirac spinors and polarization vectors. Then the basis vector labels can be 
simply {s}m = {sa, Sh,si, . . . , Sm} and the basis vectors can be 

^{{s}m)xM^,...,x^ = ^t^k^s^ • • • ■ (7.3) 

Then we have an orthonormal basis: (|{s'}m|{'S}m) is 1 if the spin labels are all the same 
and zero otherwise. 

We use a treatment of color that is conventional but more subtle than the treatment 
for spin. We turn to this subject in this section. 

7.1 Color basis 

We first note that as far as color is concerned, an initial state quark is equivalent to a final 
state antiquark and an initial state antiquark is equivalent to a final state quark. Thus, 
in the prose description in this section we use "quark" and "g" to refer to a final state 
quark or an initial state antiquark and we use "antiquark" and "g" to refer to a final state 
antiquark or an initial state quark. 

We next note that the amplitude |M({p, /}„)) us always invariant under an overall 
rotation of all of the parton colors. Thus what we really need is a basis for the space of 
color singlet amplitudes in the color space. There is a widely used and intuitively appealing 
way to do this that, furthermore, matches with the idea of color strings forming between 
outgoing partons [27] . 

The color basis vectors |{c}m,)(, are labeled by a color string configurations {c}m- A 
color string configuration can be described as a set {5*1, . . . , 5'„,} of one or more strings S. 
There are two types of strings, open strings and closed strings. An open string is an ordered 
set of parton indices that we denote hy S = [li,l2, ■ ■ ■ , ln-i,ln]- Here h is the label of a 
quark. In is the label of an antiquark, and h, ■ ■ ■ , In-i are labels of gluons. A closed string 
is an ordered set of at least two parton indices that we denote by S = (/i, • • • i ^n-i, 
Here all of the indices label gluons and we treat sets that differ by a cyclic permutation of 
the indices as being the same. Thus a complete color string configuration for a quark, an 
antiquark, and five gluons might be {[1, 6, 3, 7], (4, 2, 5)}. This is a notation for a possible 

^''Recall from Eq. (3.2) that in the case of an initial state parton, fi and /; denote the opposite of the 
incoming flavors of the parton. 
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Now we can define the basis states. We take ^{{c}m) to be a product 

^({c}„)"i'-'"™ = ^(5i)^">w ^(52)^''^P] . . . . (7.4) 

Here we have denoted the set of color indices represented in string k by 

{a}[k] = {ai,,...,aij (7.5) 
if string k is [h, ... ,ln] or {h, ... , In). 



Figure 6: Illustration an open string color basis state (left) and a closed sting color basis state 
(right). 

We can now define the component factors V{S). For notational convenience, we sup- 
pose that li = i so that the partons along the string are numbered sequentially, [1, 2, . . . , n] 
or (l,2,...,n). 

We first consider an open string, as illustrated on the left in Fig. 6. We define 

M/(5)i"} = n(5)-i/2 [f^H^^ . . . , (7.6) 

where the t"' are the SU(3) generator matrices for the fundamental representation and we 
take the ai, matrix element of the matrix product of the generator matrices (normalized 
to = 1/2). The normalization factor n{S) is 

n{S) = NcC^-^ . (7.7) 

With this normalization, 

(5|5> = ^|^(5)^'^>P = 1 . (7.8) 

{a} 

For a closed string with the same parton labels (now all gluons) we define, as illustrated 
on the right in Fig. 6, 

y^f^S)^"^ =n{S)-^/^Ti[t''H^^ ■■■t""] , (7.9) 

where, again, the are the SU(3) generator matrices for the fundamental representation^^ 
and where 

n{S) = C^ . (7.10) 



^^One could use the adjoint representation here. However, an adjoint representation string is approx- 
imately equivalent to two fundamental representation strings and having two strings when one would do 
makes the description more complicated. 
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With this normahzation, 



{S\S)^Y.MS)^'''\' = l-(^T~' . (7.11) 



{a} 



This is approximately 1 in the hmit of a large number of colors. 

A general color basis state is a product of string states. It includes a normalization 
factor n({c}m)~^/^, 

n{{c}m) = n{Si)n{S2)---n{SK) • (7.12) 
The normalization of the states is 

{{cU\{c}m) = ll{Sk\Sk) , (7.13) 

k 

where the factors are given in Eqs. (7.8) and (7.11). Thus ({c}m|{c}m) ~ 1 in the large 
Nc limit. 

The basis vectors |{c}m) are not exactly normalized and they are not orthogonal. 
However, the inner product between two different basis vectors is small in the hmit of a 
large number of colors, 

{{c'}m.\{cU) = 0{1/N^) {c'}m + {C}^ . (7.14) 

For instance, suppose {c}m = [1,2,3,4] and suppose {c'}m = [1,3,2,4], with the posi- 
tions of the two gluons reversed. Then ({c'}m|{c}m,) = — l/(-^c ~ 1)- The calculation is 
illustrated in Fig. 7. 





Figure 7: Inner products for color basis states. The left hand picture illustrates the inner product 
of the [1,2,3,4] with itself. The color diagram shown gives CpN^, which is just canceled by the 
normalization factor l/n{S) from Eq. (7.7). The right hand picture illustrates the inner product 
of the [1,2,3,4] with [1,3,2,4]. The color diagram gives — Cf/2. Multiplying by l/n{S) gives 
-l/iN^-l). 



^^The building blocks used are the invariant matrices Sij to connect a quark and an antiquark and t^j 
to connect a quark, antiquark and gluon. One could also use tlie completely antisymmetric matrix eijt 
to connect three quarks or three antiquarks. However, we don't need states made using eijt because the 
amplitudes A'l have net baryon number zero. 
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7.2 Parton insertion operators 

It will prove useful to define certain operators that act on an m parton color state with 
partons with labels {a, b, 1, . . . ,m} and add one parton with label m + 1. We let a\^{l) 
insert a giuon just to the right of parton I on whatever string contains parton /, 

\{---,[ai,---,l,---,an],---}) = \{. . . ,[ai, . . . ,l,m + 1, . . . ,an], . . .}) , 
a^+{l) |{. . . , (ai,. . ,a„), ...}) = |{. .., (ai, ... ,^,m + 1, . . .a„),. . .}) , 

(7.15) 

a\{l) \{...,[l,...,an],...}) = \{---[l,m + l,...,an],...}) , 

al{l)\{...,[ai,...J],...}) = . 

The first equation applies to the case that I labels a gluon on an open string while the 
second equation applies to the case that / labels a gluon on an closed string. The third 
equation applies when / labels a quark at the end of an open string. The fourth equation 
applies when I labels an antiquark at the end of an open string. In this case, there is no 
place to put the gluon, so the result is defined to be zero. Similarly, we define an operator 
aL(Z) that inserts a gluon just to the left of parton I on whatever string contains parton I, 



(7.16) 



We define an operator aq{l) that breaks a string at the position of a gluon with label I, 
creating a quark that we take to inherit the label I and an antiquark with the new label 
m + 1. If Z is a gluon on an open string, this creates two open strings. If Z is a gluon 
on a closed string, this turns the closed string into an open string. If I labels a quark or 
antiquark, we define the result to be zero. Thus 



alil)\{...,[ai,. 


■ ■ ,1, ■ ■ ■ , dn], ■ 


■•}> = 


K- 


. , [ai, . . 


.,m+ l,/,...a„],. 


•■}> 


aL(0|{.--,(ai,. 


. . ,1, . . . , On), . 


■•}> = 


K- 


.,(ai,. 


. .,m+ 1,1,... an),. 


•■}> 


al{l)\{.. 


.,[/,..., Un] , . 


•■}> = 


, 








aUl)\{-. 


• , K, • • • ,z], . 


••}> = 




• , [ai, ■ ■ 


/],...}> . 





aJ(Z) {...,[ai,. 


..,1,.. 


■ ) 0,n] 1 ■ ■ 


••}> 


= K-- 


.,[ai,...,m + l],[l,...an\, 


...» 


4(0 |{- ■•,(«!, •• 




■ ) 0,n)i • • 


••}> 


= K-- 


.,[/,... a„,ai,...,m + l],. 


••}> 




.,[1,.. 


■ ) 0,n] , ■ ■ 


••}> 


= , 






4(0 1{- ■ 


.,[ai,. 


..,/],., 


••}> 


= . 







(7.17) 



Finally, we define an operator al{l) that removes a gluon with label I from its string, and 
creates a new open string consisting of just a quark that we take to inherit the label I and 
an antiquark with the new label m + 1. If Z labels a quark or antiquark, we define the result 
to be zero. Thus 

4(0 !{•••> [«i> ■ ■ ■ ,1, ■ ■ ■ , an], ■ ■ ■}) = \{. . . ,[ai, . . . , an\,[l,m + 1] ■ ■ ■) , 

4(0 K- • • 1 (ai> ■ ■ ■ ,1, ■ ■ ■ ,an), . . ■}) = \{. . . ,{ai, . . . ,an),[l,m + 1], . . .}) , 

4(/)|{...,[/,...,a„],...}> = , 

4(0 |{..., [ai,..., /],... }> = . 

We will see the utility of these operators presently, when we study the color flow in parton 
splitting. 
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7.3 Color evolution for the quantum states 

In this subsection, examine how the color vector changes with successive parton splittings. 
Consider starting with a color state I?/') for m partons and letting a final state parton with 
label I and flavor fi emit a gluon with label m + 1.^^ After the gluon emission, parton / 
retains its label / and its flavor. This produces a new color state 



We can define the color splitting operator tj{fi 
component notation 



(7.19) 

fi + g) precisely by writing this out in 



aa,ab,ai,...,a; 



,...,am,am+i — 'y ^ j^^y^^°"i+l^aa,ab;Qlv 



.1; vi'Jm 



(7.20) 



The matrix T{fi)""^+^ here is the generator matrix in the 8 representation if fi = g, in the 
3 representation if fi G {u, d, . . . }, and in the 3 representation if fi G {u,d, . . . }.^^ 

This notation applies for an initial state splitting too since, from Eq. (3.2), in the 
case of an initial state parton, // and fi denote the opposite of the incoming flavors of the 
parton. In this case, the arrow in g and more generally in //—>/; + fm+i refers 

to backward evolution. 

Now consider what happens when the splitting operator is applied to one of the color 
basis vectors, |{c},„). For the case that = g we have 



m+l, a; 



(7.21) 



0.1 



Thus t\[g ^ g + g)|{c}m) is a normalization factor times a difference of two basis vectors 



n{{c+]m+i) 
ii{{c}m) 



1/2 



\{c+} 



m+l , 



n({c-}rn+l) 

n({c}„) 



1/2 



m+l , 



(7.22) 



In {c-i-}m+i, gluon I is replaced by gluons / and m + l. Specifically, if in {c}m gluon I 
appears in an open string as [ii, . then in {c+}m+i this string is replaced by 
[zi, m + 1, i„] and in {c_}m,+i this string is replaced by [ii, m + 1, i^]. 
The normalization factor is 

n({c±}m,+i) 



n({c}r 



(7.23) 



Analogous remarks apply for a closed string. We can state this result using the operators 
defined in Eqs. (7.15) and (7.16), 



(7.24) 



This is illustrated in Fig. 8. 

■^^In [14], the final state labels are called i and j. We identify i with I, and j with m+l. The labels for the 
initial m partons are {a, b, 1, . . . , m} and for the subsequent state after splitting are {a, b, 1, . . . , m, m + 1}. 
The partons not involved in the splitting keep their labels. 

^^Specifically, T(g)"j = T{q)1j = —f^i and T{g)'^^ = iftac, where q here represents any quark flavor 
and q represents any antiquark flavor. 
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Figure 8: Illustration of Eq. (7.24). A gluon I spits into gluon I and gluon m + 1. The result is 
to place the new gluon m + 1 either to the right of gluon /, with insertion operator a^(Z) or to the 
left, with insertion operator a^_{l). The normalization factors v^f ^rc not indicated. 

For a quark splitting to a quark plus a gluon we have a similar operator, one more 
generator matrix is inserted at the quark end of the string, so that 



t]{q^q + g) = y/C^a\{l) . (7.25) 

For an antiquark splitting to an antiquark plus a gluon, one more generator matrix is 
inserted at the antiquark end of the string, so that 

tl{q^q + g) = -^/C^aUl) . (7.26) 

In these equations, q represents any quark flavor and q represents any antiquark flavor. 

The operators for gluon emission obey an important identity. The matrix element M 
and our approximations to it are always color singlets, 

^ ^ J'|'^^^°"i+i^aa,ab,ai,---.Qiv.am _Q (7 27) 

I di 

Thus 

E4(//^// + g)|V'> = . (7.28) 
/ 

We will use this identity in the following sections to rearrange the splitting formulas. 

The last possibility for parton splitting is that of a gluon splitting into a quark- 
antiquark pair with flavors q and q. The color splitting operator is defined by 

IV') =4(g^<7 + g)|V') , (7.29) 

where 

^as,,a^,,ai,...,am,a,m+l = ^ ^ ^'ia,ab,ai viS; ,...,0™, (7 30) 

Our convention here is that the quark inherits the gluon label, I, and the antiquark gets 
the new label, m + 1. To analyze this, we use 



iai,am+i ^c,d ~ 2 ^c,am+i^ai,d ^jy ^"'h<^m+i^c,d ■ (7-31) 
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Thus ti{g — > g + (7) |{c}m) is a combination of two basis vectors, 
"n({ci}m+i)l I > 1 \n{{c2}, -^"1^/^ 



n{{c}m) 



m+l/ 



2N, 



im+l) 



n{{c}r 



{C2U+1) ■ (7.32) 



In {ci}m+i, the string in which gluon / resides is spht. Specifically, if in {c}m, gluon I 
appears in an open string as [ii, . . . ,iA,l,iB ■ ■ ■ in], then in {cij-m+i this string is replaced 
by two strings, [zi,...,z^,m + l][l,iB,--- ,in]- On the other hand, if in {c}m, gluon I 
appears in a closed string as (n, . . . , i^, Z, is, . . . in), then in {ci}m+i this string is replaced 
by an open string, [/, z^, . . . , i^, . . . , + 1]. In either case, 

ra({ci}m+i) _ iVc 
niicU) ~ ■ ^ ^ ^ 

In the second color state in Eq. (7.32), |{c2}m+i)i gluon / is simply removed from its 
string and a new trivial qq string is created. Specifically, if in {c}m gluon / appears in an 
open string as [ii, . . . ,iA,l,iB ■ ■ - in], then in {c2}m+i this string is replaced by two strings, 
[ii, . . . , iA,iB, • ■ ■ ,in] and [/, m + l]. On the other hand, if in {c}m gluon I appears in an 
closed string as {ii, . . . ,iA,l,iB, ■ ■ ■ in), then in {c2}m+i this string is replaced by a closed 
string {ii, . . . ,iA,iB, ■ ■ ■ in) and the open string [/, m + l]. In either case, 

n({c2}m+i) Nc 



n{{c}m) Cf 

We can state this result using the operators defined in Eqs. (7.17) and (7.18), 



(7.34) 



We can now see the advantage of this choice of color basis, beyond the fact that it is 
widely used for NLO calculations. In this basis, the description of parton splitting is very 
simple. For splittings q ^ qg and q — > gg, we simply add a gluon to a string. For splittings 
g gg, we replace one gluon on a string by two. There are two terms corresponding to 
the two possible orders for the two gluons on the string. For a splitting g — > qq, we split a 
string or we remove the gluon from its string and create a new qq string. The normalization 
factors may be considered to be just a matter of bookkeeping. However, one should note 
that the splittings q ^ qg, q ^ qg- and g — > gg come with numerical factors y/Up in 
the amplitude. These factors are large for large N^. The splitting g ^ qq comes with a 
numerical factor for the first term, [Nd {4:C-p)Y^'^ that is not large in the large limit. 
(The numerical factor for the second term is small in the large limit.) Thus the color 
factor for the splitting g ^ is smaller than the color factors for the other splittings. 
This makes g ^ qq splitting somewhat disfavored even though there are several qq flavor 
choices available. 

8. Evolution for the statistical states 

We have seen how parton splitting works at the quantum amplitude level. We now need 
to use these results to formulate the effect of splitting on the density operator p and thus 
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the statistical state |p). Suppose that we have the function 

{{P,f,s',c,s,c}m\p) (8.1) 

describing m final state partons plus the two initial state partons. After one splitting, we 
will have m + 1 final state partons, described by 

{{P,f,s',c,s,c}m+i\p) ■ (8.2) 

We need this function within the soft and collinear splitting approximations that we have 
used for the amplitudes. 

Any of the partons can split. In particular, parton i in |M({p, f}m)) and parton j in 
(^M{{p, f}„i)\ can split. In the simplest case (and the only case incorporated into typical 
parton shower Monte Carlo event generators), i = j. However, one can have i ^ j and still 
get the same momenta and flavors {p, / j^+i from both splittings and still get a logarithmic 
divergence if one were to integrate over {p}m+i- This happens when parton i emits a gluon 
m + 1 in |M({p, /}m)) and parton j emits gluon m + 1 in (^M{{p, f}m)\- An interference 
graph of this sort has no collinear divergence from Pm+i being collinear to pi or pj. However 
it does have a soft divergence from Pm+i 0. For this reason, we need to include the case 

For each splitting, we relate {p, f}m+i to the starting momenta and flavors {p, f}m 
and splitting variables {Cpj Cf} using one of the mappings 

{pJ}m+l=Rl{{pJU,{Cp,Ci}) , (8.3) 

as specified in Sec. 4. In the case i ^ j, we can use either the mapping with I = i or the 
mapping with I = j?^ We find it most useful to average over these two possibilities. We 
use the mapping with I = i with a weight Aij and the mapping with / = j with weight Aji. 
In general, the weights can depend on the momenta {p}m+i- In that case, the momentum 
mapping Ri with I = i is used for Aij and the momentum mapping Ri with I = j is used 
for Aji . In this paper, the default value is Aij = 1/2, but our notation allows for other 
choices. In the case i = j, we use the mapping with I = i = j. 

We are thus led to write the density operator p after splitting as a sum of contributions 
/5-^\ where the superscript indicates the treatment of the kinematics and the subscripts 
indicate which partons split, 

p{{P, f}m+i) = Yl PihiP^ /}m+l) 
I 

+ ^{^.,(M,n+l) f^{{p. f}m+l) + M{p}m+l) ({p, f}m+l)} 

(8.4) 



I 



Aik{{p}m+l) pf,J{{p,f}m+l)+pii{{P,f}m+l) } • 



■^^Either mapping suffices because in either case the map becomes the identity map in the soft limit 

Pm + l 0. 
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Here the sums run over the set {a, b, 1, . . . , m}. We use Eq. (6.4) for the contributions pfj 



in the case i 



p[^{{pj}m+l)=t]{fi 



fl + L+l) Vi\{p, fU+l) P{{P, fU) 



xVi{{p,fU+,)ti{fi 
xSi{{fU+i) . 



fl + fm+l) 



(8.5) 



Here Si{{f}m+i) is a counting factor that is determined by the conventions we have used 
to label the final state partons, 



Sli{f}m+l) 



f 1/2, Zg{1, 
1, 
0, 

1, 
0, 

I 1, le{a,h} 



• m}, fl = fm+i = g 

.m}, fl / g,/m+l = g 
.^}, fl = g,/m+l ^ g 

fl = qJm+1 = q 
fl = qJm+1 = q 



i.6) 



This factor is zero for // = g and fm+i = q or fm+i = q, even though flavor conservation 
allows this combination inaf;— >g + gor^— >g + g splitting, because we have chosen 
to label the final daughter gluon in these cases as parton m + 1. Similarly, this factor is 
zero for fi = q and fm+i = q, even though flavor conservation allows this combination 
in a g — > (7 + g splitting, because we have chosen to label the final daughter antiquark 
in this case as parton m + 1. In the case of a final state g ^ g + g splitting on line /, 
the splitting probability is symmetric under interchange of the labels / and m + 1 , so that 
integrating over pi and Pm+i would count the same physical configuration twice. The factor 
1/2 corrects for this symmetry. This issue is discussed in some detail in Appendix B. 
For i / j, we can use the simpler splitting operator V^°^^ as in Eq. (6.33), 



p^{{p, /}„,+i) = tlif, ^f, + fm+i) F/'^°"({p, fU+i) pUp, fU) 



t,soft/ 



X Vf°'\{p, fU+l) tjifj ^ /, + fm+l) 



(8.7) 



Note that V^'^°^^ and Vj"^^ both vanish if parton m + 1 is not a gluon. Thus /5-'^ 

i ^ j unless parton m + 1 is a gluon. 

.(0 



for 



This definition for p^- defines a mapping of the statistical states in which j/j) becomes 



\Pij ) - "^ij \P) ' 



after splitting. Then 

{{p, f,s',c',s, c}m+i\Slj'\p) 



(1)1 ^ 1 



[d{p,f, s',c',s,c}r 



X {{P, f, s', c', S, c}m+l\Sl^^\{P, /, s', C, S, C}m) 
X {{P,f,s',C,S,c}m\p) ■ 



(8.9) 
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The splitting operator 



5^ • is defined by giving its matrix elements, 



{{P, f, s',c',s,c}m+i 



^ij \{Pj,s',c',S,c}m) = 

({c', c}m+l I Q{i,j; {/}m+l) I W, c}m) 

X {{§', s}m+l\W{i, j; {/, p}m+l) I {s', S}m) 

X {{pj}m+l\ri\{p,f}m) 



(8.10) 



nc{a)nc{b) rjarjh fa/AiVa-^ lJ^F)fb/B^''T^^ ^f) 

nc{a)nc{h) f]aj)\, fa/AiVa., I^F)fb/Bi'nb, fJ-p) 




The simplest case occurs for i = j with anything other than a final state g g+g splitting. 
For that case 



({«', s}m+l\yV{l, I; {f,p}m+l)\{s', S}m) 

= Sl{{f}m+l) /}^+i)|{s}^)({s'}m|Vi({^,/}^+l)|{s'}m+l> ■ 



{{s', s}m+l\yV{l, I; {f,p}m+l)\{s', S}m) 

= Sii{f}rn+l)\^{{s}rn+l\Vi\{pJ}rn+l)\{s}rn){{s'^^^ 

+ e{le {!,..., m}, fi = fm+1 = g) {{s', s}m+i\m, I; {pU+i)\{s', s}m) } . 

(8.12) 



The second term is included in the special case of a final state g ^ g + g splitting. In 
this case, the first term is symmetric under an interchange I ^ m + 1 between the two 
final state gluons. It gives a leading singularity when gluon m + 1 is soft but also gives 
a leading singularity when gluon / is soft. We seek to use the freedom to assign labels 
to ensure that gluon m + 1 can be soft but not gluon I. The operator W in the added 
term is antisymmetric under / -^-^^ m+ 1, so that it gives zero contribution after integration 
over the final state momenta. The simplest choice for W would be the first term times 
[2 6(pm+i -Q < PrQ) — M which is antisymmetric under I m+ 1. Then the total W would 
be just the first term times 2 9{pm+i ■ Q < Pi ■ Q)- Clearly this eliminates the singularity 
when gluon I is soft. However, we adopt a shghtly more subtle procedure. Write the tensor 
that defines the three gluon vertex as v = vi + V2 + v^, 



For i = J in general we define 



Vi{Pa,Ph,Pc)=g°''{Pa-PhV ) 
V2'^"'{Pa,Pb,Pc) = g^^iPb - PcT 1 
Va'^"' {Pa ,Pb,Pc)= 5^" {Pc-Paf ■ 



(8.13) 
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Then we can define partial vertex functions analogous to those in Eq. (6.31) by 

Vj,li{p, f}m+l,Sm+l, Si, Si) 

= ^ ea{Pm+l,Sm+i; Q)* E ji{pi, Sf, Q)* e" {PU Si] Q) (8-14) 

2pm+l -pi 

X Vj'^'^{pm+i,Pi,-Pm+i - Pi) D^^{pi + p„i+i;ni) . 
for J £ {1, 2, 3}. Finally we define partial operators V by the analogue of Eq. (6.6), 

{{s}m+l\Vji{{p, f}m+l)\{s}m) = { JJ ^sA Vj^l{{p, f}m+l, Sra+1, k, Si) . (8.15) 

\j(^{Lm+l} ) 

Now the first term in W comes from Y^j k=i '^J ^ '^*K- '^^^ term with a leading singularity 
when gluon m + 1 is soft is V2 ^ The term with a leading singularity when gluon / is 
soft is ?;3 X ^3 . Therefore we need to get rid of ^3 x f 3 and double f 2 x to make up for 
the factor 1/2 in 5/. We thus define 

({s', s}m+i I >V(/, /; {p}m+l)\{s', s}m) = 

\Vl{{pJ]^+l)\{s]^){{s']m\V2,l{{pJ] m+l)\ {s }m+l / 
- {{s]m+l I Vli{{p, f]m+l )\{s]m){{s]m\ Vs,/ ({p, /}m+l) | {s'}m+l ) • 

This is computationally very simple even if it takes some time to explain. 
For i ^ j, we use the simpler splitting operators V^°^^, 

{{§', s}.m+l\yV{i,j; {f,p},n+l)\{s', S}.m) = {{s}„i+l\V^'''°^\{p, f}m+l)\{s}rn) 

X {{s'U\V^^"'\{p,f}rn+l)\{s'U+l) 



i.l6) 



i.l7) 



In the next factor, the matrix element ({p, f}m+i I'PiliP, f}m) is a delta function that 
enforces the requirement that the momenta and flavors {p, f}m are related to the momenta 
and flavors after the splitting by the mapping Qi. The definition is 

^ / [d{p,f}m] {{p,f}m+i\Vi\{p,fU)h{{p,fU)=h{{p'j'U) , (8.18) 

Ul. J 

where h({p, f}m) is an arbitrary test function and 

{{p', f'U, {Cp, Cf}} = Qi{{p, fU+i) (8.19) 

as specified in Sec. 4. Another useful identity for Vi is^^ 



1 



(m + 1)! 



[d{p,f}m+l] g{{P,f}ni+l) {{P, hm+llVllip, f}m) 

1 r (8-20) 

7—Y) E dCp9{CpGmpu,C{))9{{pJ'u+i)) , 



^®To derive Eq. (8.20) from Eq. (8.18), we add one more integration to Eq. (8.18), in the form 
/ [d{p, /}m+i] g{{p, /}m+i). Then on the right hand side, we change integration variables from {p', 
to {{p, /}m, {Cpi Cf}} according to Eqs. (4.37) and (4.70). This gives a result that is equivalent to Eq. (8.20). 
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where ^({p, is an arbitrary test function and 

iP', f'U+l = Rl{{{p, f}m, {Cp, Cf}}) (8.21) 

is the inverse transformation to Qi, as specified in Sec. 4. In this form, we display an inte- 
gration over the splitting variables that would occur in an implementation of this formalism 
as a computer program. 

The counting factor (m + 1) is the ratio of the factor (m + 1)! in the normalization 
integral for (^{p, f,s',&,s,c}m+i\<S^^j^\p) to the factor ml in the normalization integral for 
({p, /, s', c', s, c}m\p) . This factor is derived in Appendix B. The factor with parton distri- 
butions comes from Eqs. (3.17) and (3.18). In the case of a final state splitting, this factor 
is 1. 

We can now assemble our result. From Eq. (8.4) we have 

\p)=S\p) , (8.22) 

where the total splitting operator is 

S = Y,Si , 5, = 5« + ^ A, {Sll^ + 4? } . (8.23) 

I k^l 

Here the sums run over the set {a, b, 1, . . . , m} and Aik is the operator on the space of 
statistical states that multiplies a basis vector /, s', c', s, c}m+i) by the corresponding 
function Aik{{p}m+i)- 

We manipulate the result a bit. Because the quantum amplitudes are color singlets, as 
reflected in Eq. (7.28), when fm+i = g the gluon emission operators Q{1, k; {/}m+i) obey 

Y^g{l,k;{fU+i)=0 , ^g{k,l;{fU+i)=0 . (8.24) 
k k 

Then the 0(1, k; {f}m+i) are not independent. For this reason, there can be color coherence 
cancellations that are always present but are not evident if we use all of the possible 
Q{1, k; {f}m+i) operators. Accordingly, when fm+i = g, we eliminate {f}m+i) with 

the replacement 

g{l, I; {/}^+i) = "2 E {^Wi) - 2 E ^' Wi) • (8.25) 

k^l k^l 

We will see in Sec. 11 how this allows cancellations to occur at a low level of a calculation. 
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With this replacement, the total splitting operator is given for the case fm+i = g by 

{{P, f, S,C, S, c}m+l\Sl I {p, f, s',C, S, C}m) 

fa/A{fla.,IJ'F)fb/B('nh,IJ-F) 

= [m + l){{p, f}^+i\Vi\{p, f}m) —— — T~— ^ 7 — 7 2T7 — 1 2T 

^ X] \{{c,c}m+l\G{l,k;{f}m+l)\{c',c}m) 

A;e{a,b,l,...,m} 

X Aiki{p}m+i){{s,s}^+i\Wil,k;{f,p}m+i)\{s',s}^) (8.26) 

- ^{{s' , s}m+l\yV{l,l; {f ,p}m+l)\{s' , S}m) 
+ {{c',c}m+l\G{kJ; {f}m+l)\{c',c}m) 

Alk{{p}m+l) {{S, S}m+l\y^{k, I; {f,p}rn+l)\{s', s}r 

- ^{{S , s}m+l\yV{l,l; {f ,p}m+l)\{s' , S}m) 

In the case that fm+i 7^ g, that is {fi,fm+i} = {q,q}, this becomes 

{{P, /, S, C , S, c}m+l\Sl\{p, /, s', C, S, C}m) = (m + 1) {{p, f}m+l\'Pl\{p, f}m) 

nc{a)nc{b) r)s,r]]^ fa/AiVa, t^'F)fh/B^''^^^ ^f) 

nc{a)nXh) fj^fp^ fa/AiVa, fJ-j?)fb/B{r]h, fJ-p) (8.27) 

X {{c,c},n+l\G{l,l; {f}m+l)\{c,c}m) 

X {{S,s}m+l\y^{l,l;{f,p}m+l)\{s',s}m) ■ 

In the following section, we will use S to define the splitting operator at shower time t, 
7ii{t). Then in Sec. 10 we will study the evolution of the color structure that is expressed 
in Eqs. (8.26) and (8.27). In Sec. 11, we will examine the how quantum coherence for soft 
gluon emission is contained in Eq. (8.26). 

9. The operator ni{t) 

It remains to define the operator 7Yi(t). The integral of this operator over t gives the total 
probability for a splitting at any scale. 



/•oo 

/ dt Hiit) = S . 
Jo 



(9.1) 



To get the probability for a splitting at the scale corresponding to shower time t, we simply 
need to insert a delta function that defines t. Our default choice is 

Ql 



t = be 



).2 



Here the virtuality is defined using pi + Pm+i for a final state splitting and pi — Pm+i for 
an initial state splitting. We take Qq to be the hardness scale of the initial hard scattering 
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that starts the parton shower, so that the initial value of t is zero. One could take Qq 
to be the minimum of the values 2pi ■ pj for final state particles i,j from the initial hard 
scattering. Here we can neglect all quark masses compared to Qq. Other definitions of t 
are possible. For instance, many authors use a measure of the transverse momentum in a 
parton splitting. 

When we use the definition (9.2) of t, we obtain the corresponding definition of TCi^t). 

{{P, /, S, C , S, C}m+1 \ni{t)\{p, /, s', C, S, C}rn) 

= ^ {{pJ,s',C,S,c}m+l\Sl\{pJ,s',C,S,c}ra) 

Zs{a,b,l,...,m} 

(9.3) 



X 6[t — log 



Ql 



where Si is defined in Eqs. (8.26) and (8.27). In these equations there are parton distribu- 
tion functions evaluated at a factorization scale fip, which we define according to Eq. (3.38) 
in terms of t. In addition, there is a factor Og, which needs to be evaluated at a scale /xr 
that has not been made explicit in the notation. The argument presented in Sec. (3.4) 
indicates that the momentum scale at the splitting is just the resolution scale that we use 
for np. Thus we take 



10. Color evolution of the statistical states 



We are now in a position to say something about the color structure of the statistical 
states and the evolution of this structure. Notice that in general ({p, /, s', c', s, c}m|p) can 
be non-zero for {c'}m / {cjm- Even if we start with {c'}„ = {c}n, at some early stage 
of evolution, splitting will generate states with {c'}m 7^ {c}m at later stages. However, 
in the end we measure something that is color independent. This means taking the color 
trace of p. That is, we multiply [{p, f, s' , c' , s, c}m\p) by ({c'}m|{c}m) and sum over the 
color configurations {c'}m and {c}m, as in Eq. (3.34). Recall that when {c'}m / {c}m, the 
inner product ({c'}m|{c}m) is of order l/N^ to some power and is thus small in the large 

limit and numerically small for A'c = 3. Thus configurations with {c'}m 7^ {c}m are 
not very important. Furthermore, if we start with a state with {c'}m 7^ {c}m, splitting 
cannot generate {c'}m+i = {c}m+i at the next stage. Thus one can speak of a "leading 
color approximation" in which contributions with {c'}m 7^ {c}m are always dropped. This 
is what happens in most parton shower Monte Carlo programs. We do not drop terms, but 
we should understand what happens in the leading color approximation. 

In order to understand splitting in the leading color configuration, suppose that we 
apply 7ii{t) to a state \{p, f, s', c', s, c}m) in which {c'} = {c}. Thus we consider the action 
of the color splitting operator G{i,j'-, {f}m+i), on a state |{c', c}m) with {c'} = {c}, 

{{c,c}m+l\G{i,j; {f}m+l)\{c,c}ni) = ^({c}m+l |4(/i ^ /i + fm+l)\{c}ni) 

m.+l \t]{fj fj + fm+l)\{c}m)* ■ 



- 59 - 



We first consider gluon emission. That is fm+i = g- Then, after using Eq. (8.25), all of the 
contributions to 7ii{t)\p) come from Q{i,j; {/}m+i) with i / j. In order to have a definite 
case in mind, let us suppose that i and j are gluons. Then, using Eq. (7.24) to represent the 
t'^ operators in terms of the operators that insert the added gluon in particular places 
in the string basis states, our matrix element becomes 

CY{{c]m+l\a^+{i)\{c]m) ({c | ( j) | {c}m)* 

+ CY{{c}m+l\a'^-{i)\{c}m) ({c'}^+i | fll (j) | {c}^>* 

- CY{{c]m+l\a\{i)\{c}m) ({c | al (j) | {c}m) 

- CF({c}m+i|al(i)|{c}m) {{c'] m+l\a\{j)\{c] m)* ■ 

Consider the case that gluon j is just to the right of gluon i along a color string. Then 
if the new gluon is inserted to the right of gluon i, we get the same state as when the new 
gluon is inserted to the left of gluon j. That is, the third term, with a)j^(i) and a^(j), gives 
a non-zero matrix element when {c'}m+i = {c}m+i- This contribution is thus kept in the 
leading color approximation. The other three contributions would be thrown away in this 
approximation. This is illustrated in Fig. 9. 




Figure 9: Illustration of Eq. (10.2). The first diagram shows the color configuration corresponding 
to the operator G{i,j', {f}m+i)- the gluon's color is emitted from parton i and absorbed on gluon j. 
Using Eq. (7.24) to write this in terms of color basis states, one contribution is shown in the middle 
diagram. This is the leading color contribution. Another contribution is shown in the right-hand 
diagram. This contribution is subleading. If i and j had not been color connected to each other, 
there would have been no contribution that survives in the leading color approximation. 

Similarly, if gluon j is just to the left of gluon i along a color string, then the con- 
tribution with al(i) and a\j^{j) makes a leading color contribution, while the other three 
contributions would be thrown away in the leading color approximation. 

Suppose now that gluon j is not next to gluon i along a color string. Then all four 
terms would be thrown away in the leading color approximation. 

We can summarize this by saying that there is a term that is kept in the leading color 
approximation when gluons i and j are "color connected" : next to each other along a color 
string. An analogous analysis leads to the same conclusion if one or both of partons i and 
j are quarks. 

One may say that in the leading color approximation, the operating units are color 
dipoles, consisting of partons that are next to each other along a single color string. Within 
this approximation, a gluon may be considered to carry the 3x3 representation of SU (3) 
instead of the 8 representation. Then, for instance, the 3 half of a gluon forms a dipole 
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with the 3 part of the neighboring gluon or with a neighboring quark. A dipole can emit 
a gluon. But there is no interference between diagrams in which the gluon is emitted by 
different dipoles. This is illustrated in Fig. 10. The color dipole picture was introduced as 
the basis of the parton shower program Ariadne [28]. 

With the formalism presented in this paper, one 
can easily implement this approximation, but one can 
also keep corrections to it. 

We have considered gluon emission. The case of a 
g — > g + g splitting is rather different. Then the con- 
tributing terms are from Q{l,l;{f}m^i). For the op- 
erators t|(g q + q) in Eq. (10.1), we use Eq. (7.35). 
In the leading term, with operator aJ(Z), the color 
string containing the gluon breaks at the position of 
the gluon, with the creation of two new string ends. 
There is a second term, with operator al{l), in which 
gluon / simply disappears from its string and a new 
q-q string is created. However, the coefficient of al{l) 
is small in the Nc ^ oo limit and would be dropped 
in the leading color approximation. In fact, the coeffi- 
cient of the leading aq{l) is finite in the Nc ^ oo limit, 
instead of growing with Nc- Thus g ^ q + q splitting 
is color suppressed compared to g ^ g -|- g splitting 
cancelled by the number of available qq flavors. Parton shower Monte Carlo programs 
normally include the leading term in g ^ q + q splitting. With the formalism presented in 
this paper, we can keep both terms in g ^ q + q splitting. 



11. Soft gluon coherence 

When a soft gluon is emitted from the partons involved in a hard scattering, the angular 
distribution of the emitted gluon reflects the structure of the emitting partons as a whole. 
In particular, a soft gluon emitted from a pair of colored partons at an angle that is 
bigger than the angle between the partons effectively sees just one parton with the total 
color charge of the pair. Since the splitting operator Tli^t) coherently sums the leading 
soft gluon singularities for emissions from a hard partonic system, including all of the 
interference diagrams, soft gluon coherence is automatically included. It is of interest to 
see how this happens and what form it takes in the color basis used in this paper. 

Consider the matrix element of the splitting operator Si for the case fm+i = g, 
as given in Eq. (8.26) with the default choice Aik = 1/2. There is a sum over in- 
dices k € {a, b, 1, . . . , m} with k ^ I. There are two terms with different color factors, 
{{c',c}m+i\Gil,k;{f}m+i)\{c',c}m) and {{c' ,c}m+i\G{k,l; {f}m+i)\{c' ,c}m) ■ These have 
analogous structures, so it suffices to analyze one of them. The coefficient of the color 




Figure 10: The starting diagram 
from Fig. 9 in the leading color 
approximation, in which gluons are 
emitted from color dipoles. The i- 
j color dipole is a unit that emits a 
gluon and absorbs it on the other side 
of the diagram. At the next stage 
there will be an i-(m + l) color dipole 
and an (m + color dipole. 

Some of the color suppression is 
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factor -{{c',c}m+i\Gil,k;{f}m+i)\{c',c}m) is 

{{s, s}m+l|>Vtot(^, t,{f,p}m+l)\{s', s}m) = 

{{§', sU+i\W{l, I; {f,p}m+i)\{s', sU) (11.1) 

- {{S , S}m+l\y^{l,k] {f ,p}m+l)\{s' , S}m) ■ 

The first term is from the square of the amplitude for emitting the gluon from parton 
/. The second term represents the interference between the emission of gluon m + 1 from 
parton / and the emission of this gluon from parton k. Let us evaluate this in the soft-gluon 
approximation, V ysoit^ (jgfingfj Eqs. (6.35) and (6.36). We can also approximate 
Pi ~ pi and ^ Pk- With these replacements, we have 

{{s',s}m.+i\y^tot{l,k; {f,p}m+i] 

-i,Sm+i;Q)*-pi 




S],Sj S-,S^ 




Pm+1 -Pl 

l,Sm+i;Q)-Pl e{pm+l,Sm+i;Q)-Pk 



(11.2) 



Pm+l-Pl Pm+l-Pk 

There are singularities when Pm+i is collinear with pi and when Pm+i coUinear with p^. The 
singularity for Pm+i collinear with is not strong enough to give a logarithmic divergence 
when we integrate over the direction of pm+i- The singularity for pm+i collinear with pi is 
strong enough to give a logarithmic divergence. Here, the leading singularity comes from 
the first term, while the second term, representing the interference graph, gives only an 
integrable collinear singularity. 

When Pm+i is not close to being collinear with pi or pk, both terms are important. 
However, suppose that the angle between pi and pk is small and that Pm+i makes an angle 
with either of them that is substantially greater than this angle. Then the two contributions 
approximately cancel each other. That is, soft gluon radiation associated with this color 
factor is approximately confined to a cone about the directions of partons / and k with 
opening angle on the order of the angle between pi and p^. 

Thus there is approximate angular ordering of soft gluon emissions. One can also make 
an exact statement about angular ordering [5] . If we sum over the spins of the soft gluon, 
we can use 

Eu/- " ^\ i//" ~ au I Prn+lQ" ^ Q'^Pm+l Pm+lPm+l 
e''{pm+l,Sm+i;Q)e [Pm+l,Sm+i;Q) = -Q^ + : y:. j-. y^^-^ ■ 

Pm+l-Q {pm+l-Q)'' 

Sm + l 

(11.3) 

Then we can integrate over the azimuthal angle <j) specifying the rotation of Pm+i about pi 
in the rest frame of Q. This gives 

27T ^ ^Sm+i,Sm+i{{s',s}m+l\Mot{l,k;{f,p}m+l)\{s',s}m) = 

(11.4) 




soft 
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where 



^soft 



J 27r Pm+l 



Pi iPm+l-pi Pm+l-Pk 



pf I Pm+lQ'^ + Q^Pm+l Q Pm+lPm+l 



(11.5) 



Pm+l-Q {Pm+l-QY 

One can multiply this out and perform the integral using 

1 1 



27r + B"^ + 2ABcos(t) \A^ - B'^\ ' 
This gives 

J Q-Pl J Pm+l-UJ-pkPl , J 

-'soft — : 7^ S „ „ + J- 



(11.6) 



Pl-Pm+lQ-Pm+l \ {Pm+l-i^-PkY + m'^{fk) {uJ-pm+lY j {PVPm+lY 

(11.7) 

Here (5 is the velocity of parton /, 

pi-Q 

We have written uj ■ v for the product oo^v^ of a vector v with the tensor oj defined by 

Q-pi 

The tensor Pf^uj^'^ is the unit antisymmetric tensor in the pi-Q plane and projects onto 
this plane. It obeys Pf'^uj'^ai^'^i, = Pu, where Plf is the projection operator onto the pi-Q 
plane. 

In the case that all of the masses are zero, Eq. (11.7) becomes 

r Q-Pl J Pm+l-UJ-Pk ,A , . ni 1n^ 

/soft = : yr-. < I, == + 1 > (masses vanish) . (11.10) 

Pl-Pm+lQ-Pm+l [y/{pm+l-u;-pk? J 

The ratio in the first term in Eq. (11.10) is either +1 or —1. It is +1 if 

Pm+i-Q ^ Pk-Q 
Pm+i-Pi Pk-Pi 

That is, the ratio is +1 if, in the rest frame of Q, Pm.+i makes a smaller angle with pi than 
does pk- In this small angle region, we have 

/soft = — — ^ ^ (masses vanish, small angle) . (11.12) 

PVPm+l Q-Pm+l 

In the wide angle region, where Pm+i makes a larger angle with pi than does pk, we have 

-^soft = (masses vanish, wide angle) . (11.13) 
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Thus, in the massless case, after a sum over spins and an average over the azimuthal angle, 
the soft radiation outside a cone centered on pi and extending out to pk cancels exactly. This 
phenomenon is known as angular ordering. It can be used to define the angular ordering 
approximation, in which one simply neglects the radiation outside of this cone even for 
a fixed set of spins and a fixed azimuthal angle. This angular ordering approximation is 
commonly used in parton shower Monte Carlo event generators. With the formalism of 
this paper, one can make this approximation if desired, but it is not required. 

In the massive case, it is no longer true that soft radiation outside a cone centered on 
pi and extending out to pk cancels after integrating over the azimuthal angle. 



12. Inclusive evolution 



As described in Sec. 3.5, the Sudakov exponent is constructed from the operator V{t). In 
turn, V{t) is constructed from (^l\7ii{t)\{p, f, s' , c' , s, c}m) , which tells the total probability 
for the state \ {p, /, s', c', s, c}m) to split at time t. In this section, we determine the structure 
of {l\ni{t)\{p, f, s' , c' , s, c}m) . 

We know the matrix elements of TCi^t), so we simply insert the completeness relation 
(3.28) for the basis states and use Eq. (3.44) for the inner product of (l| with a basis state. 
Thus 



{l\ni{t)\{p,f, s',c',.s,c}m) 



[d{p,f,s',c',S,c}m+l] 



{m + 1)! 

X {{s}m+l\{s}m+l) {{c}m+l\{c}m+l) 

X {{P, f, s', c', S, C}m+1 \ni{t)\{p, f, s', c', S, C}m) 



(12.1) 



We use Eq. (9.3) to express 7ii{t) in terms of a sum of splitting operators Si and then use 
Eqs. (8.26) and (8.27) for the matrix elements of Si. We encounter 



1 

(m + 1)! 



[d{p, /}m+l] {{P, f}m+l\Vl\{p, f}m) 



(m + 1) 



J dCp9{Cv^mp}m,Q)) 



(12.2) 



where we have used Eq. (8.20). Thus we are really integrating over the splitting variables. 
Inside the integral and the sum over the variables {p, f}m+i are determined from {p, f}m 



and the splitting variables {CfiCf} by the transformation Ri, Eqs. (4.28) and (4.63). Then 
{l\ni{t)\{p,f,s',c,s,c}m) = 

I Ct6*,(/;) 



ric 

X 



(a)nc(6) ?7a??b fa/AiVa.^ t^F)fb/B(^^^ I^f) 



(12.3) 



nc{a)nc{h) fjs,fib fa/A{r]a, fJ'F)fb/Bir]h, Mf) 
X |6'(/m,+l/g) ({c'}m|5ii({/}m+l)|{c}m)({s'}m|?A'zK{/;P}m+l)|{s}r 

+ e{fm+l =g)Yl {W}ni\gik{{f}m+l)\{c}m) 
k^l L 

X {{s}jn\Alk{{p}m+l)Wlk[{f,p}m+l) " ^ W^ii ({/> P}m+l) | {s}m) 
+ {{<^}m\gkl{{f} m+1 

X ({s'}m|Afc(Mm+l)w^fcK{/!P}m+l) " ^ «^iK{/> ^}m+l) | {s}m 

Here the color dependent function is 

{{c}m\gij{{f}m+l)\{c}ni) = ^ {{c}^^+i\{c}m+l) 

{&,c}^+i (12.4) 
X ({c , c},„+i |^(i, j; {/jm+i) |{c', c}m) 

and the spin dependent function is 

{{s']m\wij{{f,p]m+l)\{s]m)= ^ {{s ] m+l\{s] m+l) 

{s-',s}m+i (12.5) 

X ({s',s}m+l|>V(i,i;{/,p}m+i)|{s',s}m) . 

Let us look at the color factor in Eq. (12.3) first. Using the definition (8.11) of 
Q{1, k; {f}m+i), we have 



/D 



{{c}m\gij{{f}m+l)\{c}m) = ^ {{c'}m\tj{fj ^ fj + fm+l)\{c} m+l] 

{c',c}„i+i 

X ({c},„+i|{c}m+i)^({c}m+i|4(/i ^ /i + /m+l)|{c}m) • 

(12.6) 

Then using the completeness relations (3.23) and (3.24) we find 

{{c}m\gij{{f}m+l)\{c}m) = {{c}m\tj{fj ^ f j + fm+l) t\{fi ^ fi + fm+l)\{c}m) ■ (12.7) 

The operator tjt\ is written as Tj ■ Ti in the work of Catani and Seymour on the dipole 
subtraction scheme for next-to-leading order calculations [14] . In that work, there is a sum 
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over final states, while in this paper we follow the evolution of the exclusive final state 
that comes between the operators tj and tj. We get back to the inclusive case when we 
form (l l'Hi(t) I {p, f, s',c', s, c}m) ■ The operators gij have some simple properties. From the 
definition Eq. (7.20), we see that 

{{c}m\gij{{f}m.+l)\{c}ni) = {{c}rn.\gji{{f}m+l)\{c}rn) ■ (12.8) 

Furthermore, when i = j, the operators are proportional to the unit operator. 



{{c}m.\9ll{{f}m.+ l)\{c}m) = {{c}m.\{c}m) X < 



Cf , {fijm+i} = {q,g} or {^,g} 

Ca , {flJm+l} = {g,g} 

Tr , {fiJm+i} = {q,q} 

(12.9) 

The last of these cases occurs in Eq. (12.3). 

Let us look next at the spin dependent factor in Eq. (12.3). In order to do this, we 
need to introduce the possibility of averaging over the azimuthal angle of parton splitting. 
For the splitting of parton I, we define the transverse part, q± of pm+i by 

Pm+i = api + bni + q± , (12.10) 

where ni is the lightlike vector defined in Eq. (5.3) and q± ■ pi = q± ■ n-i = 0. We let <j) 
be the angle of q± as measured in any convenient coordinate system, so that J dcp means 
integrating over q± at fixed \q±\. Thus integrating over (j) is part of integrating over the 
splitting variables ^p. 

Consider first the spin dependent factors for the case i ^ j, which arises from interfer- 
ence diagrams. Using the definition (8.17) of W{i,j;{f,p}m+i) and the orthogonality of 
the spin basis vectors, this is 

{{s'U\w,j{{f,pU+l)\{s}m) = Yl ({^'}-|^/°''(te/}-+l)|{^}-+l> 

{S}m + 1 (12.11) 

x({sWi|y/'^°^*({p,/Wi)|{4m) . 

Using the definition (6.35) and (6.36) of V^"^^, this is 

{{s'}m\wij{{f,p}m+l)\{s}m) 

= 5Z ( n h.,sjs^,s'\{^7ras) 

{s}m+i \ n.e{a,b,l,...,m} / 
^ £{Pm+l,Sm+i;Q)*-Pi e{Pm+l,Sni+l]Q)-Pj 
Pm+l-Pi Pm+l-Pj 

^ ( TT e I f4 ) "S^ g(pm+l,Sm+i;Q)*-p» g(Pm+l , gm+l ; Q) "Pj 

\ rae{a,b,l,...,m} / Sm+i 

/( n \r \ \(A e{pm+l,Sm+i;Q)*-Pi e{pm+l, Sm+l] Q) 'Pj 

= \{s }m |s}m)(47ras) } . 

' Pm+l-Pi Pm+l-Pj 

(12.12) 
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Notice that {{s'}m\wij{{f , p}m+i)\{s}m) is proportional to a unit matrix ({s'}m|{s}?n) in 
the spin indices and. that it is symmetric under i <— >■ j. The coefficient of ^{s'^tt^II^s^tt^) is a 
product of sphtting functions and Vj'^^^* as defined in Eq. (6.36), summed over Si and 
Sj and averaged over Si and Sj. That is 

{{s'}rn\wij{{f,P}m+l)\{s}ra) = {{s'}rn\{s}ni) Wij{{f , p}.m+l) , (12.13) 



where Wij is 



Wij{{f,P}m.+l) = Wji{{f,p}m.+l) 

\ E E E E E Sm+1, Sj, Sj)* ^^i°^*(te /}m+l, Sm+1, Si, Si) 



4 

Sm + l ^3 



(12.14) 



Consider next the case i = j= /G{l,..., m} with fi^g or with / = a or b and any 
fi- Using the definition (8.12) of W{l,l; {f ,p}m+i) and orthogonality for the spin basis 
vectors, this is 

{{s'}m\wil{{f,p}m+l)\{s}m)= ^ '^'Kl/lm+l) {{s'}rn\Vl{{p, f}m+l)\{s}rn+l) 

X {{s}m+l\Vi\{p, f}m+l)\{s}m) ■ 



(12.15) 



Using the definition of V as given in Eq. (6.6), this is 



fm+lj 



[{s'}ni\wu{{f,p}m+l)\{s}m)= ^ I 11 ^■Sn,sJs„,s'„]Sl{{f}r 

{s}m+i \n^{l,m+l} / 

X Viiip, f}m+l, Sm+l,Sl, s'l) Vi{{p, f}m+l,Sm+l,Sl,Sl) 



n ^s'„,s)jS,{{f}r 
Si \n^{/,m+l} / 



tm.+l) 

i{l,m+l} 

X Vi{{p, /}m+l, Sm+1, S/, s'l) Vi{{p, f}m+l, Sm+1, Si, Si) . 

(12.16) 

This is proportional to the unit matrix in the spin indices Sn,s'^ for all n except possibly 
for si,s'i, the indices that appear in the functions vi. However, in Eq. (12.3) there is an 
integration over the azimuthal angle (p as part of the integration over (^p. Once we sum 
over the final state spin indices si, Sm+i and integrate over (j), the result is invariant under 
rotations about the pi axis in the rest frame of Q. Thus the result must vanish for si s\. 
In addition, the parity invariance of the splitting vertices implies that the result is invariant 
under a reflection through a plane containing pi, Q, and any vector transverse to pi and 
Q. Under this transformation, si — s; and s'^ —s[. Thus the result is proportional to 
5si^s[- The coefficient of the 5^^^/ can be obtained by setting — > s;, summing over s;, and 
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multiplying by 1/2. Thus 



— {{s'}m\wil{{f ,P}m+l)\{s}m) 

zn 



27r 



vn = l / Sm + l Si Si 



(12.17) 



n=i / Sm+l Si 

X vIUp, f}m+l,Sm+l,Sl, S/) f}m+l,Sm+l, Si, Si) 

= {{s'}m\{s}m) Sl{{f}m+l) 

^^^\M{pJ}m+l,Sm+l,Sl,Si)\^ . 

Sm + l Si Si 

In the last line, we have not written the average over cp because, once we have summed 
over all of the spins, the result is independent of (j). 

For the special case i = j = I G {1, . . . , m} with fi = g, there is an extra term in the 
definition (8.12) of W{l,l;{f ,p}m+i)- The extra term, involving W defined in Eq. (8.16), 
is built from some of the separate terms in the three gluon vertex. Their treatment is 
essentially the same as the treatment just given for the other vertex functions. The result 
for i = j = I in general is 



7^ ({s'}m|u^«({/>P}m+l)|{s}m) = ({s'}m|{s}m)'W//({/,p}m+l) , (12.18) 

Ztt 



where 



Wll{{f,p}m+l) = Si{{f}^+i) ^ E EE 

Sm + l Si Si 



X |p«({P,/}m+l,Sm+l,S«,S/)| 

+ 9{l G {l,...,m},// = fm+i =g) 

X \v2,l{{p, f}m+l,Sm+l,SuSl)\'^ - \v^^i{{p, f}m+l, Sm+l, h, Sl){' 

(12.19) 

Here V2^i and v^ i are defined in Eq. (8.14). 

We conclude that (l|'Hi(t)|{p, /, s', c', s, c}m) has the form given in Eq. (3.49), 

(l|Wi(t)|{p,/,s',c',s,c}„) =2({s'}„|{s}„>({c'}™|/i(t,{p,/}^)|{c}„> . (12.20) 

It is proportional to the unit matrix in spin but is not proportional to the unit matrix, or 
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even diagonal, in color. The matrix {{c'}m\h{t, {p, f}m)\{c}m.) is 

{{c}ni\h{t,{p,f}m)\{c}m) = 



2 

x<5(t-log' ^0 



m + i-lY'-^+^^'^Pm+l)^ - m^fi)\ 

nc{a)nc{b) rjs,'!]-^ fa/AiVa., I^F)f'b/B^^^^ I^f) 

nc{a)nc{b) f^i^fl^, fa/AiV^, fJ'F)fb/B{llh, Mf) 
X \ G{frn+1 / g) ({c }m|{c}m) Wu{{f,p} 

m+1 ) 



(12.21) 



m \gik{{f}m+l)\{c} 1 



X [2Aik{{p}m+l)wik{{f,p}m+l) -Wu{{f,p}m+l)] 

Here {{c'}m\gij{{f}m+i)\{c}m) is given in Eq. (12.7) while Wij{{f ,p}m+i) for i / j is 
given in Eq. (12.14) and wii{{f ,p}m+i) is given in Eq. (12.19). The scales fip used in the 
parton distribution functions and ;Ur used in as are given by Eq. (9.4). Inside the integral 
and the sum over I, the variables {p, /}m+i are determined from {p, f}m and the splitting 
variables {CpjCf} by the transformation Ri, Eqs. (4.28) and (4.63). 

Notice that the matrix {{c'}m.\h{t, {p, f}m)\{c}m) is not diagonal in color. However, 
the matrix elements with {c'}m 7^ {c'}m are suppressed by powers of l/Nc- For the term 
with a factor ({c'}m|{c}m)) this was already noted in Eq. (7.14). To see this for the term 
involving gik, we should write the matrix element of gik in the form of Eq. (12.6), 

{{c'}m\gik{{f}m+l)\{c}m) = ^ ^({c j^+i |4(/fc ^ A + /m+l)|{c'}m) 

{c',c}m + l 

X {{c'}m+l\{c}m+l) jj{{c}m+l\tl{fl ^ fl + /m+l)|{c}m) . 

(12.22) 

The leading contribution for 1/Nc comes when {c'}m = {c}m and partons / and k are 
color connected: they lie next to each other on a string. As discussed in Sec. 10, we can 
use the representations (7.24), (7.25), (7.26), and (7.35) for the operators tj and to see 
that the action of tj and t[, can produce {c'}m+i = {c}m+i- If W}m 7^ {c}m, then the 
action of tj and t^, always produces {c'}m+i / {c}m+i, so that we get a color suppressed 
inner product {{c'}m+i\{c}m+i)- 



13. End of the shower 

As the shower progresses toward smaller and smaller resolution scales fi^ = Qq^"*, there 
must come a point at which the perturbative basis of the evolution equation is no longer 
valid. Then, at some evolution time tf , the shower evolution should be stopped. In the event 
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that the resolution scale of the desired measurement function is larger than fif = Qq^^*' , 
whatever happens beyond that is not seen by the measurement. Then one could simply 
apply the measurement function, calculating 

{F\p{tf)) . (13.1) 

Let us suppose that the measurement function does not see the spins or colors of the final 
state partons. Then, as in Eq. (3.34), we need 

(F|p(tf)) = / [d{pJ,s',c',s,cU]F{{pJU) {{s'U\{sU){{c'U\{c}m) 

m ' (13.2) 

X {{P,f,s',c',S,c}m\p{ti)) . 

It is significant that, although the possibility of {c'}m / {c}m is included in the shower 
evolution presented in this paper, the matrix ({c'}m|{c}m) is almost diagonal, with off- 
diagonal matrix elements being suppressed by factors of The situations with respect 
to spins and colors are different. In the end, we must have {s'}m = {s}m exactly. However, 
the spins get shuffled at each stage of shower evolution and there is no reason that a 
state with {s'}m 7^ {s}m a-t an earlier stage of evolution cannot evolve into a state with 
{s'}m = {s}m at the end. On the other hand, color differences between {d}m and {c}m 
are a little like entropy. Once {c'}m 7^ {c}m at an early stage of evolution, we can never 
get {c'}m = {c}m at the end. For this reason, the most important part of color evolution 
is the part that maintains {c'}m = {c}m throughout. This is evolution in the leading color 
dipole approximation. 

Now, what if we wish to use a measurement function with a resolution scale smaller 
than nf. Then we need a model for what happens at smaller resolution scales (or later 
and earlier proper times than given by x'^ ~ Our model should certainly include 

hadronization. We can easily extend the formalism presented here to encompass hadroniza- 
tion. We have only to replace (F|p(tf)) by 

(Fh|^/(cx),tf)|p(tf)) . (13.3) 

Here U{oo, tf) represents a model for what happens after Monte Carlo time tf . It starts with 
partonic states and maps them into the space of hadronic states. Typically the hadronic 
states are labeled by momenta and hadronic flavors but not spins. Then (Fh| represents 
the measurement function in the space of hadronic states. For purposes of discussing the 
partonic shower, we can denote 

(Fh|W(oo,tf) = (FeffI . (13.4) 

Thus Fcs is the true hadronic measurement function translated back to the partonic level. 

Assuming that the hadronization model does not use color or spin information, the 
measured cross section then takes the form 

(Fh|Z^(00,if)|p(tf)) = (Feff|/?(tf)) 

= E^ / Kp./>s'.c',s,c}m]Feff({p,/},„) (13.5) 

X {{s'}m\{s}rn){{c}m\{c}m) {{P, f, s' , C , S, C}m\p{tf)) ■ 
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Typically, hadronization models do use color information. The color field interacting 
with the outgoing partons is represented as a classical color string. The string then frag- 
ments into hadrons. This applies directly in Pythia [2] and in a different way in Herwig 
[1], where the color strings fragment into color singlet parton clusters immediately. The 
formalism of this paper is set up with the color string picture in mind. The color states 
{c}m exactly map onto string configurations, as explained in Sec. 7.1. Thus for terms in 
Eq. (13.5) with {c'}m = {c} m, one can use {c}m, as the input to the hadronization model. 
The string model does not tell us what to do with {c'}m {c}m- A reasonable suggestion 
would be to use {c'}m half the time and {c}m half the time. This could be represented as 

{Fh\U{oo,t{)\p{t{)) = Y^ I [d{p,f,s',c,s,c}„,] 

X i [Feff ({p, /, C]m) + Fefr({p, /, c]m) ] ^^^"^^ 
X ({■5'}m|{s}m)({c'}m|{c}m)({p,/,s',c',S,c}m|p(tf)) . 

If the hadronization model is based on strings, one needs to do something with the 
string ends that connect to the initial state partons. Consider, for example, the case that 
the hard collision at scale /if is a quark-quark collision and concentrate on one of the initial 
state quarks. From the point of view of backwards evolution from the scale Qq, this initial 
state quark appears as a color 3 line. The initial state quark is part of a colorless hadron, 
but the net color 3 is carried by the spectator quarks from this hadron. Thus in a color 
string model, the color strings from the hard interactions should connect to the spectator 
partons. Of course, this is not a completely simple problem. In a realistic model, the 
spectator partons from the two hadrons have many interactions with each other, possibly 
followed by their own (not very hard) parton showers. After all of these interactions, the 
spectator partons left over from a quark-quark collision must have 3 (8) 3 color, but its 
internal color state can be quite complicated. Modeling the spectator interactions is well 
beyond the scope of this paper but is addressed by Sjostrand and Skands in Ref. [12] . 

14. Conclusions 

We have presented a formulation of parton showering for hadron-hadron collisions. The 
prediction for a cross section corresponding to an observable F is given by (yF\p{tf)) or 
one of the other formulas in Sec. 13, depending on the treatment of hadronization, which 
is not covered in this paper. The dynamics of the quantum density \p{t)^ is given by the 
evolution operator Z^(t, t'), so that |p(tf)) = Z//(tf, 0)|p(0)), where |/5(0)) is determined from 
the hard matrix element that starts the shower. Thus the shower dynamics represented 
in U{t,t') is based on factorization of soft and coUinear singularities from hard scattering. 
The basic formula is Eq. (3.56), which we can rewrite as 

U{ti, t') = ATiti, t') + f'dT U{ti, t) [Hiir) - Vs(r)] AA(r, t') . (14.1) 
Jt' 

Here 7ii{t) is a parton splitting operator, as defined in the preceding sections and Vs(i) 
represents a virtual interaction that interchanges colors. The operator AA(r, t') generates 
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the standard sort of Sudakov exponential that gives the probabihty not to have an inter- 
action between shower times t' and r. We provide Table 2 to indicate where the various 
functions needed to compute U can be found. 
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Table 2: Key formulas. 



Eq. (14.1) has the proper form to conveniently generate a parton shower. Starting with 
a state at time t' , one would use M{t, t') to determine the time r for a parton splitting (or 
color rearrangement). Possibly there is no splitting before the cutoff time tf, as represented 
by the first term. Otherwise, the operators 'Hi(r) and Vs(r) give a new partonic state at 
time r. Now we operate with U{ti,T), which is to say that we apply this procedure again. 
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We offer liere concluding remarks under several headings. 

Implementation. Equation (14.1), together with the formulas in Sec. 13 and the defini- 
tions given throughout this paper, represents a certain approximation for the cross section 
cr[F] corresponding to a given observable F. When Eq. (14.1) is iterated, the result is 
expressed in the form of certain integrals and sums. It will be a significant challenge to 
find ways to implement Eq. (14.1) in a manner that allows an efficient calculation of cr[F]. 
We expect that there is more than one way to attack this problem. The choice affects the 
efficiency of calculation, but not the result, cr[F]. We leave implementation issues to later 
work. 

Evolution variable. We have chosen the evolution variable t to be proportional to 
the virtuality in the splitting. An alternative would be the transverse momentum in the 
splitting. The Herwig choice of the splitting angle does not work well with the formalism 
presented here since a parton emitted at a fixed angle can be arbitrarily soft, necessitating 
introducing the final hardness cutoff fif at each splitting. 

Momentum mapping. It is not kinematically possible for an on-shell parton to split 
into two on-shell partons. However, it is useful to approximate the mother parton as being 
exactly on-shell in calculating the (relatively) hard scattering in which the mother parton 
participates. To make this approximation, we need to take the needed momentum from 
somewhere else. Sometimes this is done by taking momentum from the mother's sister in 
the previous splitting. In the Catani-Seymour subtraction scheme for doing next-to-leading 
order calculations, there is a "spectator parton" that donates the required momentum [14]. 
Rather than taking momentum from a single parton that might not have much to give, we 
have chosen to take a little momentum from each final state parton, with each donating 
according to how much momentum it has. Note that we keep momentum in balance at 
each step, rather than waiting until the end of the shower to make adjustments. 

Common evolution. In Eq. (14.1), each parton has a chance to split or interchange 
colors between shower times t and t + dt. In the very simplest form of a parton shower, 
each parton could evolve independently, at least if one ignores momentum conservation 
and adjusts the momenta only at the end. Then the complete evolution of each parton / 
could be traced out without keeping track of the other partons. Effectively, there could be a 
separate time ti for each parton. However, independent evolution could still be implemented 
using a common evolution variable t, giving each parton its chance to split between times 
t and t + dt. The physical distinction that characterizes independent evolution is that the 
evolution is independent if the various functions involved in the splitting of parton I do 
not involve the states of the other partons. With the definitions of the operators used 
in Eq. (14.1), the use of a common shower time variable is required, since each splitting 
changes the whole partonic state and affects the probabilities for other partons to split at 
later shower times. 

Other choices. Within Eq. (14.1), there are quite a number of other choices required. 
For instance, the splitting functions must have a particular form in the limits of soft and 



- 73 - 



coUinear splittings. However, away from these limits there is a certain freedom to choose. 
Where choices like this were needed, we exercised the freedom to choose based mostly on 
conceptual simplicity. Other authors might choose differently. 

Interference and angular ordering. In Eq. (14.1), quantum interference between 
emissions of a gluon from different partons is treated exactly in the soft gluon limit. Sup- 
pression of wide angle emissions is a result. Typically, parton shower generators make 
an "angular ordering" approximation to this result. Once one has an implementation of 
Eq. (14.1), it will be interesting to make the standard angular ordering approximation and 
see how good an approximation it is. 

Spin. Partons produced in a hard scattering carry spin, which can affect the angular 
distribution of their subsequent splittings. Eq. (14.1) includes the full spin information. 
Typically, parton shower generators average over spins, thus discarding this information. 
Once one has an implementation of Eq. (14.1), it will be interesting to insert spin averages 
everywhere and see how good an approximation it is. 

Color. Partons produced in a hard scattering carry color, which can affect the pattern of 
future splittings. The formalism presented here includes the full color information. Typi- 
cally, parton shower generators make use of a leading color approximation that amounts to 
taking the first term in an expansion about l/N^ = 0. Once one has an implementation of 
Eq. (14.1), it will be interesting to make the leading color approximation everywhere and 
see how good an approximation it is. 

Improvements needed. We leave for future work the question of how one could match 
the parton shower to the exact matrix elements for 2 — > n scattering instead of simply 
starting with 2 — > 2 scattering. We also leave for future work the question of how one 
could do this at next-to-leading order. (See, however, Ref. [19].) More ambitiously, we 
would like to extend the whole formalism, including the splitting functions, to next-to- 
leading order. Current work by others on soft-coUinear effective theory may be helpful 
here [29]. 

Masses. We gave included quark masses in our momentum mappings and splitting func- 
tions. However, there are a number of issues associated with masses that we do not address. 
Suppose that we start with a hard scattering at a scale Qq that is much larger than the 
mass of the bottom quark, m-b.^'' Suppose additionally that we want to continue evolution 
down to a scale ^ that is less than m\y. Then we need a suitable variable flavor number 
scheme. At the leading order used in this paper, this is easy enough, but at higher orders 
of perturbation theory there are some subtle issues. Collins [30] has addressed some of 
these issues as they arise in deeply inelastic scattering. One should, however, note that 
for hadron-hadron collisions, the power suppressed terms that are omitted when the cross 
section is written in a factored form as a hard scattering function convoluted with parton 

^'''Here we assume that the bottom quark is the heaviest quark counted as initial state parton. One could 
include production of top quarks, but would not include the top quark as a constituent of the proton unless 
Qo > nit. 
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distributions are not of order Aqqj-,/(5o but rather of order ui'^/Qq [31]. This is not a 
problem since, in the applications we have in mind, rn^/Qo *C 1. We lack a theorem to 
tell us what to do if we want to keep the hard scattering but lower the resolution scale to 
fj? < so as to examine the final state in more detail. 

Foundations. The formalism presented here is based on ideas of factorization, both at 
the amplitude level and at the cross section level, where summation over partonic states 
that are unresolved at a given scale is essential. As the discussion of masses makes clear, 
more work is needed to make these ideas sufficiently precise to justify the formalism. 
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A. Limit on momentum fraction after splitting 

In Sec. (4.4) we defined the kinematics of initial state splitting. For a coUinear splitting 
from an initial state parton, say parton "a,", we have fja_ > rjg^. It is not exactly evident 
that this holds in away from the collinear limit, especially with masses. Here we show that 
this holds under the kinematic conditions (3.6) and (3.9). 
Let 



^(^^a) = (Pa + Ph) 

VaVhS + m (/a) + m (/b) H ; . 



^(^a) = (Pa + Ph 



i2 



= riar]hS + m (/a) + m (/b) H . 

Note that 5'(7?a) is an increasing function of its argument in the allowed region of momentum 
fractions: 

VbS (ir)a \^ fia-qhS J 

since fjaVhS > m'^ in our allowed kinematic region according to Eq. (3.9). Since ^(T^a) is an 
increasing function, we just need to show that ^(r}a) > g{f]a)- 
Recall that the kinematics requires 



(A.3) 
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where K"^ = K"^. Thus 



= {k+ pm+if -K^+ g{r)^) - g{r]^) 
= 2K ■ pm+i + m^{fm+i) 



(A.4) 



+ {m\U)-m\f,) 



1 + 



We need to show that the right hand side of Eq. (A.4) is positive. If m{fa)'^ > m{fa)'^, 
this is evident. There is only one case in which m^(/a) < m?{fa)- That is when we have a 
q ^ qg or q ^ qg sphtting in which the quark (or antiquark) enters the final state and the 
gluon enters the hard scattering. Then m?{fm+i) = "^^(/a) and m?{fa) = 0. In that case, 



(A.5) 



> 2K -pm+i - m^(/a) , 



where we have used rja^lhS > n^-'^ifh) from Eq. (3.9). The minimum value oi K -pm+i occurs 
when pm+i is proportional to K. Then K ■ pm+i = 'm{fa)VK'^, so 



This is positive as long as the condition (3.6) holds. 
We have seen that 

givn) > givs) ■ 



(A.6) 



(A.7) 



Since g{fia) is an increasing function of 17a, this implies that r)a > % 



B. Counting factors for the density matrix 

In this appendix, we organize the singular contributions to the density matrix starting with 
the quantum amplitudes defined to be symmetric in the labels of the final state partons 
(or antisymmetric in the case of identical fermions). Then we introduce the relabelings 
that define our labeling scheme for parton splittings. This produces the counting factors 
Si{{f}m+i) defined in Eq. (8.6). The counting factors are related to our parton labeling 
choices, which in turn related to the singularities of the amplitude in the limit in which 
masses can be neglected. The counting factors do not depend on parton masses. In order to 
keep our notation simple, in this appendix we simply take all the parton masses to vanish. 

Let |M({p, /}m+i)) be the exact tree level matrix element for a final state of m + 1 
partons, defined to be symmetric under the interchange of the labels for any two of the 
final state partons, or antisymmetric if the two partons are identical fermions. The matrix 
element may have a singularity when any of the dot products of two parton momenta. 
Pi -pj, approaches zero. Let us define approximate matrix elements \M{{p, f}ra+i'-,i, j)) 
that approximate the complete matrix element when Pi-pj, approaches zero. Here i G {a, b} 
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and j G {1, . . . , m + 1} or z, j S {1, . . . , m + 1} with i < j. There is more than one way to 
do this. We can, for instance, use the method of Sec. 6. The approximate matrix elements 
|M({p, f}m+i', hi)) thus defined should include a theta function 

\Vi-Vj\ < , min \i>k-k\ (B.l) 

SO that |M({j3, i, j)) is not singular when • for some other pair of partons 
approaches zero. Given the approximate matrix elements, the complete matrix element 
can be written as 

m+l m+1 

|M({p, /}„+!)> J]|M({p,/|^ +i;a, j)) + m+i;b,j)> 
i=i i=i 

■m+i (B.2) 

+ E |Af({p,/}™+i;i,j)> . 

The right hand side of Eq. (B.2) approximates the complete matrix element in any of the 
singular limits. For our purposes, it is convenient to define 

|M({p,/Wi;*,i)> = |Af({p,/Wi;i,^)> (b.3) 

for z, J G {1, . . . , m + 1}, i > J. Then we can symmetrize the third term in Eq. (B.2), 

m m 

|M({p,/Wi)) ~ Yl \M{{p,fU+i;a,j)) + Y\M{{pJU+r,h,j)) 
i=i i=i 

^ m+l (B.4) 

+ 2 E \M{{p,f}m+i;i,j)) . 
Now we construct the density operator from 

P0{{pJ}rn+l) = \M{{p,fU+l)){M{{pJ}rn+l)\ . (B.5) 

Imagine expanding both |M({p, and (^M{{p, f}m+i)\ according to Eq. (B.4). 
There are a number of terms, with the general form 

const. x\M{{pJ}^+i;i,j)){M{{p,fU+i;t',j')\ ■ (B.6) 

Not all of these contributions have soft or collinear singularities strong enough to produce 
a logarithmic divergence if one were to integrate over the momenta {p}m+i- 

One contribution that does have a strong enough singularity comes when £ 
{1, . . . , m + 1} and i' = i, f = j or i' = j , j' = i. The sum of these contributions is 



m+l 

2 \Mi{p,fU+i;i,j)){M{{pJ}m+i;i,j)\ , (B.7) 



«J=1 
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where we have used Eq. (B.3). 

Another contribution that has a strong enough singularity comes when i = i' = a and 
G {1, . . . , m+1} with j' = j. Similarly, we can have i = i' = h and G {1, . . . , m+1} 

with j' = j. The sum of these contributions is 

m+1 

J2 \M{{p, /Wi; a,j)){M{{p, /}„,+! ; a,j)\ 

, (B.8) 

m+1 ^ ' 

+ ^ \M{{pJU+i;h,j)){M{{pJU+i-h,j)\ . 
j=i 

Another contribution that does have a strong enough singularity comes when i, j, i' ,j' £ 
{1, . . . , m+1} and i' = i, j' ^ j or else j' = j, i' ^ i or else i' = j, j' / i or else j' = i, i' ^ j. 
These four cases are really the same, with different labeling. Adding the contributions, we 
have 

m+1 m+1 

^ ^ \M{{pJ}rn.+i;i,k)){M{{p,fU+i;j,k)\ . (B.9) 

i,j=l k=l 

Here we have a leading singularity only if = g. This is the interference between gluon 
emission from line i and gluon emission from line j. 

We can also have interference between gluon emission from an initial state line and 
gluon emission from a final state line or between gluon emission from one of the initial 
state lines and gluon emission from the other. When we add all of these cases together and 
add them to the contribution in Eq. (B.9), we get 

m+1 

J2 \M{{pJ}m+i;i,k)){M{{p,fU+i;j,k)\ . (B.IO) 

i je{a,b,...,m+l} k=l 

There are no more combinations of i, j, i', j' that give leading singular contributions to 
p. Thus the sum of the leading singular contributions is obtained by adding the contribu- 
tions (B.7), (B.8), and (B.IO). 

This formulation is fine for constructing an inverse shower, starting from a state with 
many partons and combining partons to reach a hard scattering with fewer partons. In 
order to construct a shower starting from the hard scattering, it is convenient to adopt 
a labeling convention in which /3({p, /j^+i) is not symmetric under interchanges of the 
parton labels. This is easy to do. 

Consider the contribution (B.7). Using the i <-> j symmetry of the matrix elements, 
this is 

m 

\M{{V. f}m+i;i, m + 1)>(M {{p, f}m+i; i,m + l)\ 

i=l 

m m (B.ll) 

+ 9 EE ^Wi; i,j)){M{{p, /Wi; . 



2 

i=i j=i 
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We can now deliberately break the relabeling symmetry in the second term by interchanging 
the labels j and m + 1. That is, we choose to label the daughter parton that here caries 
the label j by m + 1 instead. Note that this changes p{{p, f}m+i) for any fixed value 
of {Pi/lm+i- However, the result of integrating p{{p, f}m+i) against any measurement 
function (which must be symmetric under label interchanges) stays the same. After the 
interchange j <-> m + 1 , the sum over j in the second term simply becomes a factor m — 1 . 
This gives 



m 

E 

i=l 



M{{p, fU+i; i,m + l)){M{{p, f}m+i;i, m + 1)\ 



m 



^ m 



i=l 



(B.12) 



\M{{p, f}m+i]i, m + l)>(M({p, i, m + 1) | 



This is 



m + 1 



J2 \M{{p, /Wi; hm + l)){M{{p, f}m+i;i, m + 1)|; 



(B.13) 



We now consider the possibilities for flavors: {/j, fm+i} could be {g, g}, {g, q}, {g, q}, {q, g}, 
{q, g}, {q, q}, and {q, q}, where q stands for a quark flavor, q stands for an antiquark flavor, 
and a q together with a q stands for a quark flavor and its corresponding antiquark flavor. 
Thus we can insert a factor 

l=Oi{fi,fm+l}={g,g}) 

+ 0{{fi, L+i} = {g, q}) + 0{{{k fm+i} = {g, q}) 
+ e{{fi,L+l} = {q,g})+e{{fi,U+l} = {q,g}) 

+ Oiifi, L+i} = {q, q}) + 0{{k L+i} = {q, q}) ■ 

In the coefficients of the {g, q} and {g, q} theta functions, we can further define the labeling 
by interchanging the labels i and m + 1, so that m + 1 is the label for the gluon. In the 
{q, q} term, we can further define the labeling by interchanging the labels i and m + 1, so 
that m + 1 is the label for the antiquark. With these label choices, our contribution is 



(B.14) 



{m + l)J2 Si{{f}m+i) \M{{p, /Wi; l,m + l)){M{{p, fU+i; l,m + l)\ 
where 



(B.15) 



1=1 



Sl{{f}m.+l) 



1/2 
1 , 
0, 
1 , 
, 



/e{i, 



^m}, fi = fm+l = g 

fl / gJm+l = g 
.^}, fl = gJm+l + g 

. W> h = 9,/m+i = q 
h = 9, /m+1 = q 



(B.16) 



We thus derive the factor (m + 1) and the factor 1/2 for a final state g — > g + g splitting. 

Consider now the contribution (B.8). We can break the relabeling symmetry by inter- 
changing the labels j and m + 1. Then there are m + 1 equal terms, giving 

(m + 1) ^ 5K{/Wi) |M({p,/}„+i;/,m+l)>(M({p,/}„+i;^,"i + l)| , (b.17) 

ie{a,b} 
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where 

Si{{f}m+i) = l /G{a,b} . (B.18) 

Consider, finally, the interference diagrams, Eq. (B.IO). We separate this into several 
terms according to the values of i and j, 

m+l 

E \M{{pJU+i;t,k)){M{{pJU+v,j,k)\ 



ij6{a,b} fc=l 



m m+l 



+ E EE l^^({^^'^Wi;^,A:))(M({|5,/Wi;j,A;)| 



ie|a,bj j=l k=l 
m m+l 



+ E E E l^^(i^'-^Wi;i,^)>(M(te/}m+i;i,^)| 

i=l ie{a,b} fc=l 
m 

+ E El^^({^^'^Wi;^,^)>(M({p,/}m+i;m + l,A:)| (B.19) 

ie{a,b} k=l 

m 

+ E EI^^(i^^'-^Wi;"i + i,A:)>(M({p,/Wi;j,^)| 

je{a,b} fc=l 



+ E E /Wi; m + l, fc))(M({p, /}m+i; j, A:) I 
j=l k=l 

m m 

+ E E ^' k)){M{{p, /}m+i; m + 1, A;)| 



i=l fc=l 
m m+l 

+ E E |M({p,/}m+i;^,^))(^mp,/}m+i;j,^)| ■ 

i:i = l fc = l 

We relabel the indices, treating each term separately. In each case, we interchange k <-> 
m + l. In the first term, this gives m+l equal terms from the sum over k. In the second 
and third terms, this gives m equal terms from the sum over k. In the fourth through 
seventh terms, each term in the sum over k remains as one term. Finally, in the eighth 
term, there are m — 1 equal terms from the sum over k. After relabeling, we have 

(m + l) \M{{p,f}m+i;i,m + l)){M{{pJ}„,+r,j,m + l)\ 

ije{a,b} 
m 

+ m Y El^^(^^'^~Wi;^,"i + i)>(^'^({p,/}m+i;i,m + i)| 

je{a,b} j=l 
m 

+ "^E E \M{{p,f}m+i;i,m + l)){M{{pJ}m+i;j,m + l)\ 

i=l je{a,b} 
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+ Yl ^|^({p,/Wi;i,"^ + i)>(^({p,/Wi;^,"^ + i)| (B.20) 

ie{a,b} k=l 
m 

+ Yl J^|M({;),/}„^+l;fe,m + l)>(M({p,/}^+l;J•m + l)| 

je{a,b} k=l 
m m 

+ ^"Wi; k,m + l)>(Af({p, /}„+i;i,m + 1)| 

j=i k=i 

m m 

+ E E + i)>(^({P' ^, + 1) I 

i=l A:=l 

m 

+ {m-l)Y, \M{{p, f}^+i;i,m + l)){M{{p, f}m+i;j,m + 1)\ . 

These terms can be combined, after changing the names of some of the summation 
indices, and added to the contributions (B.15) and (B.17) to give the revised density 
operator, which we can call pi, 

Pi{{pJ}m+i) = 

(m + 1) Yl Si{{fU+i) \M{{pJU+i;l,m + l)){M{{p,f}rn+i;l,m + l)\ 

i£{a,b,l,...,m} 
m 

+ (m+l) \M{{pJ}„^+i;l,m + l)){M{{p,f}m+i;k,m + l)\ . 

/,fc£{a,b,l,...,m} 

(B.21) 

The first line here contains the direct terms, Eqs. (B.15) and (B.17), while the second line 
is the interference terms, Eq. (B.20). 
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